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44 Sommaire. — L'auteur expose les résultats d'expériences qui mettent en évidence l'influence de la 


_ Mécanisme de la relaxation nucléaire. — 
+ Les phénomènes d'absorption liés à la résonance 
nucléaire dépendent de l'existence d’une pola- 
…risation partielle d’un ensemble de moments magné- 
… tiques nucléaires. Ces phénomènes se produisent 
…_ par l'intermédiaire des mécanismes qui conduisent 
«à l'établissement de l'équilibre thermique, c’est-à- 
- dire qui conduisent à une répartition de population 
… entre les sous-niveaux Zeeman nucléaires donnée 
« par la formule de Boltzmann. 

Cet échange nécessaire d’énergie thermique entre 
… les noyaux et l’ensemble du système, suit, dans le 

- L 


» temps, une loi exponentielle e 7, loi qui fait inter- 
. venir la grandeur caractéristique 7, appelée « temps 
de relaxation spin-réseau », et dont les valeurs 
“peuvent varier entre 10°‘ et 10's suivant les subs- 
. tances et les conditions. 

- La figure 1 montre l'effet précédent (temps de 
… relaxation dans un cristal de NH,Br à la tempé- 
- rature de l'air liquide, T; — 325). 

… Les mécanismes de relaxation font intervenir 
… essentiellement des champs alternatifs d’origine 
- ihermique, capables de provoquer des « réorienta- 
- tions » du spin nucléaire. Dans le cas des moments 
… magnétiques purs (des protons, par exemple), ces 
… champs doivent être des champs magnétiques 


n (:) Conférence faite à la Société française de Physique 
» Je 15 janvier 1954, traduite par M. J. Brossel, École normale. 


pression sur les temps de relaxation de la résonance nucléaire des protons dans les liquides. Des mesures 
ont été faites jusqu’à plusieurs dizaines de milliers d’atmosphères. Les résultats obtenus ne peuvent 
pas s’interpréter par la seule variation du coefficient de viscosité et du coefficient de diffusion en fonction 
de la pression. Ces dernières grandeurs dépendent essentiellement des degrés de liberté de translation 
2 des molécules dans le liquide, tandis que les temps de relaxation dépendent aussi fortement des mou- 
vements de rotation des molécules. Une augmentation de pression inhibe davantage les mouvements 
de translation que les mouvements de rotation; ce qui explique que les temps de relaxation dépendent 
d’une manière différente de la pression que la viscosité et la diffusion. 


internes ayant des composantes oscillantes dans la 
zone des radiofréquences. 
Dans le cas de noyaux possédant un moment 
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Fig. r. — Mesure du temps de relaxation de NH, Br à 88° K 
par une méthode directe. 


électrique quadrupôle Q, des gradients de champs 
électriques ayant des composantes alternatives 
peuvent provoquer les transitions, mais ce dernier 
cas n’est pas envisagé dans ce qui suit. 


Origine des champs fluctuants. — C’est dans 
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les liquides que la relaxation nucléaire a été étudiée 
tout d’abord. 

Les champs magnétiques internes peuvent avoir 
deux origines 


a. Moments magnétiques des autres noyaux de 
la même molécule; 
b. Moments des noyaux des autres molécules. 


L'origine de la fluctuation est le mouvement 
thermique des molécules : le champ d’un dipole 
dépend de l'angle et de la distance. 4 

La fluctuation, dans le cas a (d’origine intra- 
moléculaire), provient de la rotation de la molécule; 
dans le cas b (d’origine intermoléculaire), elle pro- 
vient de l’arrivée ou du départ des molécules voi- 
sines. 

La distinction est importante et doit être conservée 
à l'esprit. 

La figure », montrant un ensemble de molécules 


Fig. 2. 


d’eau, permet d'illustrer simplement les remarques 
précédentes. 


Spectre dufchamp interne. — [La grandeur de 
ces champs internes est d'ordinaire de quelques 
gauss. Leur efficacité à produire les transitions 


dépend cependant de la rapidité des fluctuations 


incohérentes qui déterminent la répartition d’inten- 
sité du ‘spectre. 

Dans un liquide ordinaire, le temps caractéristique 
de réorientation ou de changement de position 
d’une molécule est de l’ordre de 10‘! s, ce temps 7°, 
le temps de « corrélation », est très court par rapport 
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à la période ®, de la précession nucléaire. Le spectre | 
est « blanc » au voisinage de ,. Mais son intensité 
en ce point est d’autant plus faible que le temps T. 
est plus court (fig. 3). 


Intensité spectrale 


w, DIRE 
Fig. 3. 


Donc, plus les fluctuations sont lentes, plus elles 
sont efficaces pour provoquer la relaxation nucléaire, 
et l’on peut s'attendre à ce que 


Relations entre la relaxation nucléaire et la 
viscosité et la mobilité moléculaire. — La 
mobilité moléculaire traduit sa présence par l’exis- 
tence de la viscosité. Des expériences dans lesquelles 
on faisait varier la viscosité dans de grandes limites 
en faisant varier la température (dans la glycérine, 
par exemple) ont montré que n est proportionnel 


x 


à T,, et par suite que T, est proportionnel à 


il 
Tr 

Il existe une vérification encore plus directe, dans 
le cas des liquides polaires, où la « dispersion diélec- 
trique » de Debye permet d'atteindre directement le 
temps de corrélation de l'orientation moléculaire. 
L'accord entre les deux types de mesures est bon 
dans ce cas. 

Il est une quantité physique qui est encore plus 
directement reliée à T,, en particulier en ce qui 
concerne le champ intermoléculaire, c’est la cons- 
tante de diffusion D. La plupart des théories de 
l’état liquide (Frenkel, Eyring, par exemple) pré- 
voient une relation entre et 1), dont la forme 
générale est semblable à celle de la relation de 
Stokes-Einstein 

AT 


6714 


D = 


Nous avons donc trois quantités dont il est 
intéressant de mesurer les variations : 


19 Temps de relaxation T, (dépendant du temps 
de corrélation 7. des champs internes); 


=. nr 7 
29 Viscosité, 2) 
l 
J - 6 : l 
30 Constante de diffusion, D' 


But des expériences à haute pression. — 


il 


2 ÿ RAS à 7 


- La plupart des liquides présentent un accroissement 
considérable de viscosité en fonction de la pression, 
par des facteurs de 20 à 5o à 10 000 kg/cm? 


x 


… (Bridgman). La raison est, à coup sûr, la réduction 


très différent de celui qui résulte d’un abaissement 
de température. 

Que deviennent la « mobilité moléculaire » et 
« l’orientabilité moléculaire » telles qu’elles nous 


M Ces questions sont à l’origine du travail que nous 
… décrivons ci-dessous (travail effectué à Harvard 
ki par G. Benedek). 

+ Le problème était de mesurer T, en fonction de P 
ct comparer les résultats à ceux de Bridgman 
relatifs à n. 

M) Au cours des mesures, il apparut qu'on pouvait 
MA aussi mesurer D par la résonance nucléaire, et qu’on 
_ pouvait par suite comparer les variations de # 
É et de D en fonction de la pression. 


_ Méthode utilisée. — La résonance nucléaire a 
été observée (et 7, mesuré) par la méthode de 
«l'écho de spin » décrite par Hahn. 

Considérons les noyaux comme un ensemble de 


30°? 180° 
! 8 


+ moments dans le plan équatorial (impulsion dite 
‘de 9o°). 

_ A-partir de ce moment-là, les moments ont une 
- précession libre dans le plan équatorial. Mais ils 
perdent leur cohérence de phase parce que le 
— champ H, est toujours légèrement différent d’un 

point à l’autre. 

… On applique alors, à l'instant { une impulsion 
« de 1809 » (c’est-à-dire de longueur double) dont le 
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« toupies magnétiques » dans un champ fort H;. 
Les axes sont animés d’une précession autour 
de H,. 

Un champ magnétique H,, perpendiculaire à H,, 


: J Amplificateur Amplificateur 
à large bande | là bandeétraite 


k 14 Ka 
"aile ‘ “echo” 
NE l'es 
PAT pulse second pulse 
Fig. 4. — Appareil pour écho de spin. 


(Pour la théorie générale, 
voir _E. L. HAN, Phys. Rev., 1950, 80, 580.) 


tournant autour de H, avec la vitesse de précession, 


fait basculer les moments magnétiques. Ce champ 
de radiofréquence est appliqué sous forme d’une 
impulsion de durée telle qu’elle amène juste les 


10 millisec. 


t 
10  millisec. 


Fig. 5. — L’homogénéité du champ magnétique est diminuée dans le 
cas de la figure du bas. 


résultat est de tourner face pour face le plan équa- 
torial de la figure, si bien qu’à l'instant 2 { la cohé- 
rence de phase, entre les divers moments, est entie- 
rement rétablie, ce qui provoque l'émission d’un 
« écho ». 

La figure 4 représente l’appareillage utilisé et la 
figure à montre une séquence d’écho. 

Les champs internes incohérents provoquent une 


baisse de l’intensité de l’écho, c’est-à-dire une baisse 
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de la cohérence de phase qu’il n’est pas possible versal » T,, mais, on peut aussi mesurer Ti par 
de rattraper par « l'impulsion de 1800 ». Ceci est, une variante de la méthode et, dans les liquides, 
en fait, un effet du « temps de relaxation trans- T, et T7, sont identiques, de toute manière. 


1 î 
10 20 30 40 50  millisec. 


f 
10° 20 millisec 


Fig. 6 (haut). — Courbe T, : Expositions multiples superposées. 
Combinaisons : go°-180°. Distance des échos : 4 ms. 

(bas). — Courbe T,: Expositions multiples superposées. 

Combinaisons : 1802-90. Distance des échos : 2 ms. Glycérine. 


10 p— : a nie 2 La partie « haute pression » de l'expérience 

comportait comme problème essentiel la construction 
d’une enceinte non magnétique tenant la pression 
et permettant l'application d’un champ magnétique 
intense. Un alliage Cu-Be fut utilisé sur la suggestion 
de Bridgman. Un cylindre de 5 cm, avec un diamètre 
interne de rem, tient 10000 atm sans subir de 
distorsion permanente. 

Un avantage de la méthode de l’écho de spin est 
qu'il suffit de faire une seule amenée coaxiale dans 
le volume haute pression, et qu’il n’est pas nécessaire 
d’avoir de champ modulateur. L'échantillon et la 
bobine haute fréquence étaient scellés à l’intérieur 
d’un silphon flexible immergé dans le liquide compres- 
seur. / 

Les techniques standard de Bridgman ont été 
utilisées, les échantillons étaient des liquides qui 
ne se solidifiaient pas à haute pression et étudiés 
par Bridgman en ce qui concerne les variations du 
coefficient en fonction de la pression. 


n pentane (CsH») d 
———n hexane (C4 Hu) \ 
—---- lodure d'éthyle (C.Hsl)f\ 
si... Joluêne (Cr Hg) 


Résultats. — En général, T; décroît quand P 
croît, ce qui indique l’augmentation attendue de %%, 


700 x DL ) 12000 ou encore, la diminution de la mobilité moléculaire. 
rene dE Cependant, dans tous les cas envisagés T, décroît 
Fig. 7. beaucoup moins que % n’augmente ,le produit Z,n 


n'est pas du tout constant mais croît avec la pression 

La figure 6 représente l'effet d’affaiblissement La figure 7 permet de comparer à cet égard Je 

de l’écho en fonction du temps £ séparant les deux comportement de quatre liquides organiques : 
impulsions de radiofréquence. pentane, hexane, iodure d’éthyle et toluène. 


— Dans l’hexane, par exemple, croît par un fac- 
teur 4o quand 7, décroît par un facteur 7 seulement. 

L'eau ne présente qu'un changement modéré de %, 
mais, là encore, 7, change encore moins. 

En résumé, il apparaît que le temps de corré- 
- lation des champs internes, ne change pas, en 
… fonction de la pression, dans la même proportion 
… que la viscosité du liquide. 


Mesure de la constante de diffusion D). — 
Dans un champ inhomogène, il y a quelquefois une 
perte de cohérence de phase dans la précession des 
spins à cause de leur migration lente d’un point 
à un autre par suite d’une diffusion. 

Ceci produit une décroissance de l'écho, d’un 
type tout à fait caractéristique, avec une loi en e * 
au lieu de e-’. La théorie complète de cet effet a 
été faite, et on peut l’utiliser pour mesurer quanti- 
…. tativement la diflusion dans les liquides, du moins 

dans quelques liquides. 

La figure 8 montre une décroissance d’écho 
(dans CH,1) caractéristique de l'effet de diffusion. 

Toutes les méthodes de mesure de D nécessitent 
un étiquetage des molécules. Dans le cas présent, 
l'étiquette qui permet de reconnaître la molécule, 
est la direction, ou la phase, du moment magné- 
tique attaché à la molécule. 

La méthode cependant possède des limitations 

: sérieuses dans ses applications, à cause du fait 
. suivant : si la molécule contient plusieurs noyaux 
semblables dans des positions non équivalentes 


Fig. 8. — Courbe de décroissance due principalement 
à l’effet de diffusion. 


(Iodure de méthyle.) 


Expositions multiples superposées : 900-1800. 


tn à 
(7) = 1,6.10°/st, D — 68 cm°/jour. 


\ 


» (comme dans l’éthane, par exemple) la courbe- 
enveloppe de l'écho est déformée par un autre 
effet résultant du couplage des spins nucléaires, 
et l'effet de diffusion est masqué. 

Dans l’iodure d’éthyle tous les H sont équi- 
valents. 
Un autre avantage de la méthode est qu’elle 
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permet de détecter aussi la convection qui affecte 
l'enveloppe de l’écho de façon très différente. 

Les résultats des mesures de diffusion sont les 
suivants : 

Dans le cas de l’eau, les changements de D et n 


—- = 


L ! 
Jodure de méthyle (CHa1) 


Valeurs relatives du temps 
de relaxation 


(7, (P)/T, ()) 


Valeurs relatives 
du coefficient 
de diffusion 


2000 40 6000 


00 
Pression (kg/cm°) 


Fig. 9. 


en fonction de P sont en bon accord. C'est-à-dire 
que Dn est pratiquement constant; comme le 
prévoit la loi de Stokes-Einstein. 

Dans le cas de CH, I, il n’existe pas de mesure 
de n. Cependant la variation de D peut être comparée 
à celle de 7. La figure 9 montre que le changement 
de D est beaucoup plus grand que le changement 
de T, ce qui montre encore que T. change moins 
vite que la mobilité moléculaire. 


Conclusion. — Une explication générale de 
ces résultats est possible, semble-t-il, si l’on se 
souvient que deux effets provoquent les fluctuations 
incohérentes des champs internes : 


a. la rotation d’une molécule; 


b. la mobilité, et, à un moindre degré, la rotation 
de ses voisines. 


Maintenant, si nous réduisons, par ‘accroissement 
de pression l’espace vide entre molécules, en entas- 
sant les molécules dans un espace moindre, il est 
raisonnable de penser que l’on réduit ainsi beaucoup 
plus les mouvements d’une molécule d’un endroit 
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à l’autre que les mouvements de rotation propre 
d'une molécule dans l'espace qui lui est alloué. 

Ceci est certainement vrai pour des molécules 
sphériques. 

La viscosité et la diffusion dépendent de la mobi- 
bilité moléculaire tandis que la relaxation nucléaire 
dépend en outre à un degré beaucoup plus grand 
des degrés de liberté caractérisant l'orientation. 

Il semble que les expériences précédentes montrent 
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cette différence. Elles montrent aussi la relation 
qui existe entre D et n. 

Elles n’en sont cependant qu'à leur début, ét 
beaucoup reste à faire. 

Enfin, notons en terminant, que la résonance 
nucléaire permet de mesurer D de façon précise 
(de façon plus précise même que les autres méthodes) 
dans les cas limites où la méthode s'applique. 


Manuscrit reçu le 20 mars 1954. 
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SUR UNE NOUVELLE MÉTHODE 
POUR LA COMPARAISON PRÉCISE DU RAYONNEMENT DES AMPOULES DE RADIUM 


Par Mme IRÈNE CURIE, 


Laboratoire Curie de l’Institut du Radium. 


Sommaire. — La méthode de compensation entre les courants produits dans deux chambres d’'ioni- 
sation se prête à des mesures de comparaison.de très grande précision entre des sources de rayonnement 
d'intensité presque égale. Par l’emploi d’une ampoule de radon, dont la décroissance est connue avec 
précision, utilisée comme étalon auxiliaire, on peut adapter cette méthode à la comparaison du rayon- 
nement d’ampoules de radium d’intensités différentes. 

On décrit des modèles de chambres propres à ces mesures, et des expériences préliminaires de compa- 


raisons d’étalons de radium., 


L'emploi d’ampoules de radon comme étalon auxiliaire a parfois de l’intérêt même dans le cas de 
mesures faites avec une seule chambre d’ionisation ou avec d’autres appareils de mesure. Il peut faciliter 
la comparaison d’ampoules d'intensité extrêmement différente, ou bien la comparaison d’étalons qui se 
trouvent dans des laboratoires différents et dont on désire éviter le transport. 


Malgré le très grand nombre de radioéléments 
artificiels qui ont été découverts, le radium en équi- 
libre avec le radon et le dépôt actif à vie courte est 
jusqu'ici le corps le plus commode pour servir 
d’étalon de rayonnement . 

Un étalon de radium est constitué par une ampoule 
contenant une quantité soigneusement pesée d’un 
sel de radium extrêmement pur. En principe, toute 
ampoule ainsi préparée est un étalon. En pratique, 
comme la préparation d’une telle ampoule est un 
travail chimique extrêmement délicat, il est impor- 
tant de comparer avec une grande précision le 
rayonnement y des tubes étalons, pour s'assurer 
que le rapport de l'intensité du rayonnement corres- 
pond au rapport des poids. D'autre part, une ampoule 
contenant un sel de radium non pesé peut être utilisée 


. comme étalon secondaire si son rayonnement a été 


comparé avec précision à celui d’une ampoule 
étalon. Pourvu qu'il n’y ait pas d’impuretés radio- 
actives, il n'est pas nécessaire alors que le radium 
soit tout à fait pur. Il est donc extrêmement impor- 
d'avoir une bonne méthode de comparaison 

des tubes de radium par le rayonnement +. 
Les ampoules utilisées comme étalon ont été 


d'abord : une ampoule préparée par Marie Curie 
en 1911, qui a été choisie comme étalon inter- 
national, et une ampoule préparée par Hônigschmidt 
à la même époque [1]. | 

Ces tubes ont été soigneusement comparés entre 
eux et le rapport des rayonnements y reconnu 
égal à 2 pour 1000 près, au rapport des poids. 
Les ampoules destinées à servir d’étalon secondaire 


+ 
4 
LE 
t 


À 


dans les différents pays étaient comparées succés- 


See à Paris à l'Étalon International, à Vienne 
à l’étalon de Hôünigschmidt, Les valeurs ont toujours 
été concordantes à 2 pour 1000 près. 


En 1934, 20 ampoules étalons ont été préparées 
par Hônigschmidt. Deux de ces ampoules ont 
remplacé respectivement l'Étalon International et 
l'étalon de Vienne : on craignait, en effet, que les 
ampoules anciennement préparées n'éclatent un 
jour, du fait de la pression de l’hélium venant du 
rayonnement «;'c’est pourquoi, du reste, on a faït 
des ampoules relativement grandes pour la quantité 
de sel qu’elles contenaient. Malheureusement, les 
ampoules de Hônigschmidt n’ont pas été comparées 
entre elles par le rayonnement +, et quelques-unes 


fl 


Mines 


eut ont été comparées aux étalons anciens. 
de Paris et de Vienne. 


qu'il était important de faire avec la meilleure 
« précision possible les comparaisons qui n'avaient 
pas été faites, d'autant plus que certains doutes 
avaient été émis sur la concordance des divers 
étalons. 

Les comparaisons d’étalons avaient été faites 
po au Laboratoire Curie de lInstitut du 
… Radium de Paris et à l’Institut für Radiumforschung 
_ de Vienne. Les appareils utilisés étaient différents, 
* mais les deux instituts employaient des chambres 
…._d'ionisation. Au Laboratoire Curie on utilise le 
… grand condensateur du service de mesures et la 
méthode de compensation au moyen du quartz 
… piézoélectrique [2] ©. 

La précision de 2 pour 1000 était obtenue en 
faisant de longues séries de mesures dont on prenait 
la moyenne. La précision des mesures nn 
est sensiblement plus faible. 

J'ai cherché une méthode qui permette de faire 
des déterminations individuellement plus précises 
afin de pouvoir comparer plus rapidement des 
_ étalons. 

On peut espérer augmenter beaucoup la précision 
de comparaisons faites au moyen des chambres 
d’ionisation en employant la méthode de compen- 
sation : Soit deux chambres d’ionisation semblables 
. portées à des tensions opposées et reliées au même 
électromètre, et des sources qui donnent dans les 
deux chambres des courants 1, et 1, peu différents; 
le courant mesuré est la différence 1 — 1, — 1. 
Si les courants 1, et 1, sont suflisamment intenses 
et réguliers, la mesure de 1, ou Z, et de it permet de 
calculer le rapport de Z, et J, avec une précision 
beaucoup meilleure que celle que l’on obtiendrait 
en mesurant séparément 1, et J,. 

On peut aussi se servir de la méthode de compen- 
+ sation pour rechercher de très petites différences de 
période : ainsi Marie Curie l’a utilisée en 192/{ pour 
essayer de mettre en évidence une différence entre 
- la période d’un tube de radon ordinaire et celle 
d’un tube de radon très concentré : elle a pu vérifier 
… que la période était la même à 0,5 pour 1000 près, 
- précision que l’on n'aurait pu atteindre par des 
mesures directes. 

Peu après, j'ai utilisé le dispositif de chambres 
. en, compensation établi par Marie Curie à cette 
. occasion, pour déterminer avec précision la période 
du radon par une méthode nouvelle : on observe 
les instants 4, &, {’ où un courant constant produit 
. par une ampoule de radium dans l’une des chambres 
- est égal aux courants variables Z,, 1,, 1! produits 
dans l’autre chambre par deux ampoules de 
radon d'intensité peu diflérente que l’on y place 


() Le modèle de la chambre est dû à Debierne et date 
de l’époque de la comparaison des premiers étalons. 


La Commission mixte de Radioactivité a estimé 
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d’abord séparément, puis ensemble. La différence 
entre la moyenne des temps # et & et le temps 
donne la période. La précision obtenue était de 0,5 
pour 1000 [3]. 

J’ai recherché si l’on ne pourrait pas trouver une 
méthode de compensation pour comparer les étalons 
de radium, mais on ne peut utiliser avec avantage 
cette méthode qu'avec des sources de rayonnement 
de valeur peu différente, ce qui n’est pas toujours 
le cas dans les comparaisons d’étalons. 

Il est possible de surmonter cette difficulté en 
comparant les ampoules de radium non l’une à 
l’autre, mais à une ampoule de radon, puisque la 
variation d'activité de celle-ci avec le temps est 
connue avec précision. 


Principe de la méthode. — Prenons deux 
chambres semblables montées en compensation, 
soit, par ordre décroissant de quantité de radium, 
des ampoules À, B, C, ., que l’on veut com- 
parer entre elles. On met dans l’une des chambres 
le tube À et dans l’autre chambre une ampoule de 
radon dont l'intensité de rayonnement est un peu 
plus forte. On mesure la différence 1, des intensités, 


um = Ti— Ir. 


x, d’abord négatif, augmente quand le courant x 
dû au radon diminue, passe par o à un temps 4, 
puis devient positif. On enlève alors l’ampoule A 
et on la remplace par l’ampoule B, plus faible. 
Le courant à, négatif d’abord, passe par o pour un 
temps {», puis devient positif, on remplace alors 
l’ampoule B par l’ampoule C, etc. 

Le rapport des intensités des tubes A, B, C, 
est égal au rapport des intensités du tube de radon 
aux temps {4, ln, Le, . La précision obtenue sur 
le rapport des intensités dépend de la précision 
obtenue sur la détermination de l'instant où le cou- 
rant s’annule : la constante radioactive du radon 
étant environ 2,1.10 6s71, on voit qu'une différence 
de 5os correspond à une différence de 104 sur 
l'intensité du courant. 

D'autre part, on peut calculer facilement que, 
la période du radon étant déterminée à 0,5 pour 1000 
près, l'erreur venant de l’incertitude sur la valeur 
de la période est égale à n fois 0,5 pour 1000, n étant 
le nombre de périodes de décroissance du radon, 
entre les temps {4, fn, Le, considérés. Si, par 
exemple, on a : 


ER Rt — 0,8 


À 


= 0.5, he 0,23, 


>| Q 
æ) 


l'erreur possible de ce fait sera de 0,15, 0,5 et 1 
pour 1000 respectivement sur ces trois rapports. 
L'incertitude sur la valeur de la période du radon 
n'intervient donc pratiquement pas dans la compa- 
raison d’ampoules d'intensité voisine et l’on voit 
que, même pour des ampoules dont l'intensité serait 


différente dans un rapport de 1 à 4, l’erreur qui 
peut provenir de là est extrêmement faible. 

Pour que la précision ne soit pas illusoire, il est 
nécessaire que la sensibilité de la chambre, par 
milligramme de radium, soit rigoureusement la 
même pour les diverses ampoules. Or les ampoules 
de radium ont parfois des dimensions sensiblement 
différentes, et nous avons vu que certaines sont 
relativement grandes pour la quantité de sel qu’elles 
contiennent (par exemple 20 mg de sel dans une 
ampoule de 3 cm de long et 3 mm de diamètre). 
Il est donc nécessaire d’utiliser des chambres peu 
sensibles à un petit déplacement de la source, dans 
le sens de la longueur et dans le sens du diamètre. 

La sensibilité par milligramme de radium risque 
aussi de varier un peu pour des sources d'intensité 
différente, si le courant d’ionisation n’est pas tout 
à fait saturé; il convient donc d’avoir des chambres 
où la collection des ions se fait le mieux possible. 
Cependant, remarquons qu'il n’est pas nécessaire 
avec la méthode ci-dessus, que la sensibilité des 
chambres soit exactement la même pour les diverses 
intensités de source : il suffit que la sensibilité 
varie avec l'intensité de la même façon dans les 
deux chambres, ce qui doit être réalisé si les deux 
chambres sont identiques. 

D'autre part, pour que l’on puisse déterminer 
avec précision les temps 4, fn, fc, ..., il faut que 
le courant i varie assez rapidement avec le temps, 
c’est-à-dire avoir un fort courant /,k dont une faible 
fraction soit aisément mesurable. Il faut donc des 
chambres donnant une valeur élevée du courant par 
millicurie, et naturellement un courant aussi régulier 
que possible. 

Il faut empêcher que la source placée pour agir 
sur une chambre ne produise un effet sur l’autre 
chambre; ceci peut être réalisé en mettant les 
chambres assez loin l’une de l’autre et en installant 
des écrans de protection convenables. 

La précision sur la détermination des temps 4, 
l», etc. ..., n'indique pas forcément la précision 
réelle des expériences. Si la sensibilité relative des 
deux chambres changeait un peu au cours des 
mesures, cela se verrait immédiatement par l'allure 
anormale’ de la variation du courant i, mais on ne 
fait des mesures que pendant une période de temps 
de l’ordre de une ou deux heures avant et après les 
temps 1, {», . . :, et si, entre ces périodes, la sensibilité 
de l’une des deux chambres changeait, ceci produirait 
seulement sur toute la série de mesures, un décalage 
d'intensité que l’on ne pourrait reconnaître. Il y a 
un moyen de vérifier la constance de la sensibilité 
relative des deux chambres : la chambre I contient 
une ampoule de radon fixe et la chambre II reçoit 
successivement les ampoules A, B, C;, ...; dans 
l'intervalle des mesures des ampoules de radium, 
on met dans la chambre II une autre ampoule de 
radon à peu près de la même intensité que celle de 
la chambre I; on mesure alors une faible différence 
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de courant qui doit décroître exponentiellement 


avec la période du radon, mais la moindre variation » 


de la sensibilité relative des deux chambres pro- 
duirait sur ce courant une variation notable, que 
lon mettrait facilement en évidence. 


Chambres d'’ionisation. — Les diverses condi- 
tions que doivent remplir les chambres ne sont pas 
faciles à réaliser simultanément. On peut toujours 
s'arranger pour que l'influence d’un petit dépla- 
cement de la source soit faible, si l’on met la source 
très loin de la chambre, mais ceci ne permet pas 
d’avoir un fort Courant par milligramme de radium. 
D'autre part, quand le courant d’ionisation est 
intense, la saturation est plus difficile à obtenir. 

La grande chambre d’ionisation de forme plate, 
dec. au service de mesure du Laboratoire Curie 
(fig. 1) est satisfaisante à divers points de vue. 


de l'ähneau de garde 


ä l'électromètre 


Fig. x. 
La boîte plate cylindrique T est reliée à la tension, le disque E 
est relié à l’électromètre. La source S est placée au milieu 
d’un disque de plomb de 1 cm d'épaisseur, de 52 em de 


diamètre, qui couvre la chambre. La distance de l’élec- 
trode isolée E à la surface supérieure de la boîte T est 
de 35 mm; cette boîte a une hauteur intérieure de 55 mm. 
L’électrode est un disque d’aluminium de 2,5 mm d’épais- 
seur, de 50 cm de diamètre, et la boîte est en laiton nickelé 
de 1,5 mm d'épaisseur. 

Le disque E est vissé sur la plaque reliée à l’électromètre : 
c’est par erreur qu’il en est isolé sur la figure, 


L'’intensité par milligramme de radium est assez 
grande, environ 0,034 u. e. s., et l’influence d’une 
différence de longueur des sources est très faible, 
car on peut déplacer une source dans un cercle 
de 7 em de rayon, sans avoir de différence de cou- 
rant supérieure à 0,5 pour 100. La sensibilité à 
une différence de position en hauteur est de l’ordre 
de 1 pour 100 pour 3 mm, de sorte qu'il peut 
être nécessaire d’en tenir compte. La saturation 
du courant est assez bien réalisée, la distance 
des électrodes étant seulement de 3,5 cm. Des 
expériences de comparaison d’ampoules seront pro- 
chainement entreprises avec deux chambres de 
Cesbype 

Remarquons qu'avec ces chambres qui ne sont 
pas étanches, il doit y avoir des variations de sensi- 
bilité au cours du temps, en raison des changements 
de température et de pression. En principe ceci 
ne doit pas avoir d'inconvénient fondamental, 
pourvu que la variation de la sensibilité soit Ja 
même dans les deux chambres; il sera donc important 
qu'elles soient bien à la même température. 


Î 


. OL 


1% Mutre part, j'ai fait construire des chambres 
selon le principe représenté dans le schéma de la 
figure 2. La source est placée au milieu d’un tube 
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Fig. 2. 


. Une enveloppe intérieure reliée au sol et pouvant supporter 
un gaz sous pression (région hachurée) contient la chambre 
d’ionisation constituée par un cylindre E, relié à l’électro- 
mètre, placé entre deux cylindres parallèles T et T’ portés 
à la tension. La distance des cylindres T et E, ou E et T'’ 
est de 3 cm. Dans le cylindre intérieur, de diamètre 5 cm, 
on glisse un tube de plexiglas P dans lequel se trouve la 
source S. 


central, entouré par une chambre cylindrique symé- 
\/ trique; on peut prévoir, qu'avec cette disposition, 
…. un déplacement de la source dans le sens de sa 
longueur ou perpendiculairement doit avoir un 
* effet nul en première approximation. La chambre 
- peut être remplie d’un gaz sous pression pour aug- 
- menter la sensibilité. 
M. Kipfer a réalisé deux chambres du modèle 
ci-dessus et publiera le dessin d'exécution ainsi 
que l’étude de l’appareillage. 
Remarquons qu’une fois les chambres montées 
en compensation, il est extrêmement facile de faire 
les corrections de déplacement d’ampoule, les essais 
de saturation, etc., en plaçant deux sources de 
-radium presque égales dans les deux chambres 
“et en faisant varier dans une seule chambre la 
position de l’une des ampoules, la tension, etc. 


…— Essais préliminaires (2). — Avec les chambres 
. cylindriques, remplies de CO, à 15 atm de pression, 
… les rayons y étant filtrés seulement par la paroi 
… de 3 mm de fer de la chambre, le courant produit 
par une source de o,1 mg Ra est environ 0,27 u. e.s. 
On utilise une tension de 1500 V; si l’on place 20 mg 
de radium dans chacune des chambres et que l’on 
ajoute ensuite la source de 0,1 mg d’un côté, celle-ci 
produit à peu près le même courant qu’en l’absence 
des ampoules, de 20 mg, ce qui montre que, pour 
. un courant de l’ordre de 54 u. e. s. dans les chambres, 
le courant est encore bien proportionnel à l'intensité 
de la source. 

On peut, sans avoir de différence supérieure à 1 
- pour 1000, déplacer une source de 3 cm de part et 
d'autre du milieu de la chambre, ou la déplacer 


(2?) Ces essais, auxquels j’ai participé, ont été faits princi- 
_ palement par M. Kipfer et par M. Lecoin, et seront l’objet 
d’une publication par la suite. 
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D dde à l’axe, de 3 mm (les ampoules 


étaient fixées dans des tubes de plexiglass placés 


selon l’axe des chambres). 


La mesure du courant compensé est évidemment 
d'autant moins précise que la compensation est 
meilleure; d’après la précision obtenue sur ce cou- 
rant pour un facteur de compensation déterminé, 
on calcule que la précision sur la mesure du courant 
d’ionisation dans les chambres est de l’ordre 
de 0,5: 10“. 

Quand on met en compensation une ampoule de 
radium et une ampoule de radon, le moment 
où le courant est nul peut être déterminé avec une 
précision d’autant plus grande que le courant dans 
les deux chambres est plus grand. La figure 3 donne 
l'exemple des mesures de courant au voisinage de 
l'intersection pour une ampoule de radium de 20 mg. 
On voit que l'intersection peut être déterminée à 
moins de os près, ce qui correspond à une erreur 
sur l'intensité inférieure à 10‘. Les mesures de 
courant étaient faites par compensation au moyen 
d’un quartz piézoélectrique. 

La vérification de la sensibilité relative des deux 
chambres a été faite en mettant une ampoule de 
radon à la place des ampoules de radium, comme il 
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F9: 

Variation-du courant à au voisinage du temps d’intersection 
qui est environ à 115 mn après le moment où l’on avait 
commencé les mesures. Le courant Z dans les chambres 
est environ égal à 54 u. e. s. On voit que l’écart des mesures 
par rapport à la droite moyenne est en général inférieur 
à 0,005 u.e.s., soit o,1 pour 1000 du courant principal. 
Au voisinage de l'intersection, la variation relative sur 
le courant à est trop rapide pour que l’on puisse faire 
les mesures, mais le moment où l’on observe le courant 
nul coïncide bien avec l’époque donnée par l'intersection 
de la droite avec l’axe d'intensité o. 

La lame de quartz donne 7,4 u. e. s. de charge électrique 
pour : kg posé sur le plateau. Le courant est exprimé 
en « grammes par seconde ». 


a été dit plus haut. Au cours d’une expérience qui 


a duré 7 jours, et au cours de laquelle 15 ampoules 


ont été comparées, la sensibilité relative des chambres 
est restée remarquablement constante. 
Les expériences d’essai ont montré que la préci- 
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sion des mesures est certainement très supérieure 
à 1 pour 1000, mais avec les étalons dont on dispose, 
la précision réelle est plus faible. Ces ampoules 
contiennent une quantité de sel de radium qui 
produit une absorption appréciable de son propre 
rayonnement ; la correction due au verre des ampoules 
n’est pas négligeable non plus. Il y a une incertitude 
importante sur ces corrections, du fait que l’on ne 
sait pas quelle est la proportion de radon occlus 
dans le sel ou réparti dans l’ampoule; la grosseur 
des grains du sel de radium a aussi une importance; 
enfin dans les ampoules de dimensions relativement 
grandes, préparées par Hôünigschmidt, le sel peut 
se déplacer beaucoup : on essaye de le répartir 
aussi uniformément que possible, mais on ne peut 
pas être sûr d’avoir toujours la même répartition. 
Si le sel est tassé à un bout du tube, l’absorption 
est beaucoup plus grande. 

Pour diminuer ces inconvénients, il faut filtrer 
le rayonnement y par une épaisseur de métal équi- 
valente à plus de 5 mm de plomb afin d'utiliser seu- 
lement les rayons les plus pénétrants du dépôt 
actif, mais même alors les corrections sont notables. 
Aïnsi pour un étalon contenant environ 24 mg de 
radium et dont la paroi est de 0,27 mm de verre, 
on a calculé une correction d'absorption du sel 
de l’ordre de pour 1000 et une correction du 
verre de l’ordre de 3 pour 1000; l’absorption du 
sel peut être encore plus forte si les cristaux sont 
relativement gros. On se rend compte qu'avec les 
étalons qui existent actuellement, même en répar- 
tissant le sel le mieux possible, on a bien une incer- 
titude de l’ordre de 30 pour 100 sur les corrections, 
et la précision réelle sur la comparaison des étalons 
ne doit pas être supérieure à 2 ou 3 pour 1000. 
Du reste la pesée du sel doit comporter aussi une 
incertitude de l’ordre de 1 pour 1000. 

Les différences entre le rapport des poids et le 
rapport des rayonnements de plusieurs tubes ont été 
au plus de 4 pour 1000, dans des expériences faites 
avec la seule filtration du tube central de la chambre. 
On pourra sans doute améliorer un peu la précision 
en utilisant une filtration plus grande. 

Un essai fait dans ce sens a donné de mauvais 
résultats, parce que la filtration -était réalisée avec 
un tube de fer épais dont le canal central était 
trop étroit. Or on peut voir géométriquement que 
si la source est trop proche du tube absorbant, la 
moindre différence dans sa position radiale change 
notablement l’épaisseur de métal traversée, de sorte 
que l’on ne peut comparer correctement des ampoules 
de diamètre un peu différent. Il est nécessaire que 
les sources soient placées dans un tube d’une matière 
peu absorbante et que le tube absorbant, pas trop 
épais, soit appliqué contre la paroi intérieure de 
la chambre d’ionisation. On pourra sans doute 
faire des mesures dans de bonnes conditions avec 
un tube de plomb de 5 à 10 mm d'épaisseur, de 
diamètre extérieur égal au canal central, soit 5 cm. 
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Gependant, le plus commode sera peut-être a 


faire des mesures au moyen de deux conden- 
sateurs plats du type de la figure 1. Naturellement » 
la précision sera moins grande, le courant étant. 
beaucoup moins intense, mais elle sera sans doute. 
suffisante étant donné la nature des étalons que 
l’on a à mesurer. 

Les nouvelles chambres cylindriques permettront 


de déterminer avec beaucoup plus de précision la 


période du radon, et elles pourront servir à comparer 
des étalons de radium avec une précision très supé- 
rieure à 1 pour 1000 quand on aura des ampoules 
permettant une telle précision. 


Possibilités offertes par l'emploi des ampoules - 


de radon dans les étalonnages. —— La méthode 


consistant à comparer entre elles les ampoules de … 


radium en utilisant une ampoule de radon comme 
étalon d'intensité variable ne permet une très 
grande précision que combinée avec la méthode 


de compensation. Cependant elle n’est pas sans 
intérêt même quand on ne dispose que d’une seule « 


chambre, sur laquelle on place alternativement 
l’ampoule de radium et l’ampoule de radon. En 
effet, on peut déterminer par des séries de mesures 
alternées les temps /,, 45, {c, où l’ampoule de radon 
donne la même intensité que les ampoules À, B, C, 


et en déduire le rapport de ces ampoules. Cette … 
détermination peut être meilleure que la détermi-” 
nation directe du rapport des intensités des courants » 


d’ionisation, parce que la méthode revient à 


comparer des courants égaux, ce qui élimine un, 


certain nombre de causes d’erreurs, en particulier 


celles qui peuvent provenir d’une saturation impar- … 


faite du courant d’ionisation, si les ampoules sont 
d'intensité assez différente. On peut aussi déter- 
miner le temps {, avec un appareil de mesure et le 


temps f# avec un autre appareil, plus sensible, si la 


source B est beaucoup plus faible que la source A. 
Ainsi, par l'intermédiaire du radon, il est possible 


1 


d’étalonner des ampoules de radium beaucoup plus ! 


faibles ou beaucoup plus fortes que les étalons 
pesés, qui contiennent ordinairement 20 ou 30 mg 
de radium. La précision de 0,5 pour 1000 sur la 


# 


période du radon permet déjà de comparer deux « 


ampoules de radium différant d’un facteur 100, 
avec une précision supérieure à 0,5 
certainement possible de déterminer la période du 
radon avec une précision de l'ordre de 1074, en 
utilisant les nouvelles chambres, ce qui augmentera 
encore les possibilités offertes par cette méthode. 


beaucoup plus fortes que les étalons, contenant par 
exemple plusieurs centaines de milligrammes de 
radium, la correction d'absorption par le sel devient 
très importante, mais bien souvent il n’est pas néces- 
saire de faire cette correction, quand on s'intéresse à 
l'intensité du rayonnement qui sort de l’ampoule 
et non à la quantité de radium qu'elle contient. 


pour 100. Il est : 


| 
Naturellement, dans le cas d’ampoules de radium 


Ÿ 


Les ampoules de radon peuvent aussi servir 
m dl'auxiliaire pour la détermination précise du rapport 
e sources de rayonnement y pénétrant autres que 
les ampoules de radium. Les conditions sont théori- 
……quement un peu moins bonnes, parce que si la 
sensibilité des deux chambres changeait un peu, 
par suite de variations de température ou pour 
… toute autre raison, la variation pourrait être un 
… peu différente si le rayonnement y est différent, 
alors qu’elle devrait être identique quand les rayon- 
US sont identiques dans les deux chambres. 
En pratique s’il s’agit de rayonnements y filtrés 
par plus de 5 mm de plomb, il est probable que 
… ceci intervient peu. 

…. Le radon peut être utilisé pour résoudre la diffi- 
.culté pratique qui consiste à comparer des étalons 
de radium qui se trouvent dans des pays différents. 
_ Le transport d’ampoules de radium est une opé- 
- ration ennuyeuse à cause de la valeur du produit 
jet des formalités de douane. Le radon, produit de 


® [1] Curie M. — J. Physique Rad., 1919, 2, 795. 
n [2] Curie M. — J. Physique Rad., 1933 4, 450. 


S nu COMPARAISON PRÉCISE DU RAYONNEMENT DES AMPOULES DE RADIUM 795 


faible valeur, peut être expédié bien plus faci- 
lement, On peut mesurer avec précision une ampoule 
de radon.par rapport à un étalon de radium, à 
Paris par exemple, puis expédier le radon par avion 
dans un autre pays, où il sera comparé à une autre 
ampoule de radium que l’on veut étalonner. 
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du Haut-Katanga, qui a assumé les frais de cons- 
truction et d'installation des chambres d’ionisation 
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physicien, M. Kipfer, de consacrer une partie de 
son temps à l’étude de ces chambres et aux expé- 
riences préliminaires de comparaisons d’étalons. 

Je remercie aussi M. Kipfer, qui a construit et 
étudié les chambres et M. Lecoïin, qui a contribué à 
établir les conditions de mesures et participé à celles- 
ci. M. Kipfer a remis à la commission mixte de radio- 
activité, en Juillet 1953, un rapport sur les premières 


comparaisons d’étalons effectuées par nous en 1953. 
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THÉORIE CINÉTIQUE DES PLASMAS HOMOGÈNES FAIBLEMENT IONISÉS. I. 


Par Micumez BAYET, 
Faculté des Sciences de Toulouse, 


JEAN-Loup DELCROIX, 
École Normale Supérieure, Laboratoire de Physique, Paris 


et JEAN-FRANÇoISs DENISSE, 
Observatoire de Paris, Section d’Astrophysique. 


— On définit une méthode rigoureuse d'intégration de l’équation de Boltzmann par 
approximations successives pour un gaz faiblement ionisé placé dans un champ électrique alternatif 
et un champ magnétique continu. On développe la fonction de distribution des vitesses électroniques 


Hohner: 


puis on développe chaque fonction f, en fonctions sphériques. 

Dans cette première partie la méthôde est appliquée au modèle de Lorentz parfait des gaz électrons- 
molécules (on néglige totalement les échanges d'énergie entre électrons et molécules). On montre que la 
méthode des fonctions sphériques réussit particulièrement bien parce que les fonctions sphériques 
sont des fonctions propres de l’opérateur intégral J qui figure au second membre de l’équation de Boltz- 
mann et qui représente l'effet des chocs. On montre que cette propriété est générale et que les valeurs 
propres s’interprètent physiquement comme étant des fréquences de relaxation des anisotropies du 
plasma, la première d’entre elles, y, coïncidant avec la fréquence des chocs classique dont la définition 
est en général assez arbitraire. On effectue ensuite le calcul des approximations d'ordre 0, 1 et 2 par 
rapport au champ électrique; l'approximation d'ordre 1 permet d'étudier le tenseur de conductivité du 
plasma; l’approximation d'ordre », ainsi que le cas du «gaz de Lorentz imparfait », seront traités dans 
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1. Introduction. — «a. NOTION DE GAZ DE 
LORENTZ. — On appelle gaz de Lorentz un gaz 
constitué de deux types de particules : des « par- 
ticules lourdes P » de masse M à raison de N par 
unité de volume, et des « particules légères p » de 
masse m en nombre beaucoüp plus restreint, soit n 
par unité de volume; par définition on a donc 


n EN, (4) 


m << M. (2) 


Un tel gaz jouit des deux propriétés fondamentales 
suivantes : 


10 par suite de la relation (1), chaque particule p 
peut être considérée comme se déplaçant indivi- 
duellement dans un gaz de particules P; 

20 par suite de la relation (2), au cours d’un 
choc p—P, il n’y a qu’un très faible échange 
d'énergie entre les deux particules. 


Nous dirons que le gaz de Lorentz est parfait si 


[1) SEAT Er 
les rapports : et sont considérés comme nuls 


ne 

N M 

de telle sorte qu'on peut négliger complètement 
l'interaction p —p et l’échange d'énergie dans un 
choc p — P. Si, au contraire, on tient compte des 


. . LL : . 
corrections du premier ordre en ;, (faible inter- 


action p —p) et en ee (faible échange d’énergie 


dans un choc p—P}), nous dirons que le gaz de 
Lorentz est imparfait. 

On rencontre des gaz de Lorentz dans l'étude des 
phénomènes physiques suivants 


— diffusion d’un gaz léger dans un gaz lourd; 

— diffusion des neutrons dans les réacteurs 
nucléaires ; 

— diffusion des électrons dans les métaux et les 
semi-conducteurs ; 

— diffusion des électrons dans les gaz ionisés. 


b. GAZ DE LORENTZ ÉLECTRONS-MOLÉCULES. 
Le présent article est consacré au cas particulier 
où les particules p sont des électrons et les par- 
ticules P des molécules à symétrie sphérique; ce 
modèle permet de représenter assez correctement 
les propriétés des plasmas de décharge : un tel 
milieu contient n électrons, n' ions positifs et 
N molécules neutres par unité de volume; on sup- 
pose d'ordinaire n — n' tout au moins en moyenne 
ce qui n’est, en général, possible que si l’on néglige 
les effets dus aux parois. En tous cas les plasmas 
qu’on étudie au laboratoire sont généralement 
peu ionisés de sorte qu’on a simultanément n — N 
et n'— N; passer de ce cas réel au modèle de Lorentz 
c'est supposer n' assez faible pour que les effets 
dus aux ions positifs soient négligeables; ceci mérite 
discussion puisque, comme nous venons de le voir, 
n et n' sont du même ordre; la discussion repose 
sur deux points : tout d’abord si nous désirons 
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étudier seulement les propriétés électromagnétiques 
du milieu, c’est-à-dire les déplacements pris par les 
charges électriques sous l'influence de champs 
électriques et magnétiques nous pouvons négliger 
les mouvements des ions qui du fait de la différence 
des masses sont d’amplitudes bien inférieures à ceux 
des électrons. Ceci dit, les mouvements des électrons 
résultent des champs appliqués et des interactions 
entre particules soient : 


— interactions électron-électron : e—e; 
— interactions électron-ion : e — P+; ] 
— interactions électron-molécule neutre : e — P, 


On est tenté de considérer comme négligeables 


les deux premiers mécanismes devant le troisième 
en s'appuyant sur le fait que l’on a 


n € N et n' ÆN; 


en fait, les particules chargées ont proportionnelle- - 


ment un effet beaucoup plus important que les 
molécules : du fait du grand rayon d’action de la 
force de Coulomb tout se passe comme si chaque 
particule chargée occupait un volume efficace bien 
plus grand que celui occupé par une molécule neutre; 
il est par suite nécessaire que n et n' soient de plu- 
sieurs ordres de grandeurs inférieurs à N pour qu’on 
puisse ne tenir compte que des interactions e — P[1]: 
Si, au contraire, l’ionisation est importante comme 
c’est le cas dans les atmosphères stellaires les inter- 
actions e—P deviennent importantes; cela ne 
veut pas dire qu’on a un gaz de Lorentz e— P* 


parce que la relation (1) n’est pas satisfaite de sorte 


que l'interaction entre les électrons est aussi impor- 
tante. 


€. LOI DE FORCE ÉLECTRON-MOLÉCULE. La 
force d'interaction e—}P sera supposée centrale; 
nous n’aurons que fort peu à préciser la loi de 
force; notons cependant que, si l’on suppose que les 


électrons polarisent les molécules proportionnelle- 


: ? x , M) ! j 
ment au champ qu'ils créent, c’est-à-dire à 3 On 


obtient une force d'interaction électron-dipole induit 
en 


une loi analogue a été considérée par Maxwell qui 
a montré comme nous le verrons plus loin qu'elle 
conduisait à des résultats particulièrement simples; 
mais la force considérée était une force de répulsion 


alors que la force électron-dipole induit est une 


force d’attraction analogue aux forces de cohésion 


I 0 s , 
en ; de Van der Waals-London qui s'exercent 


entre deux molécules [2]. Il est certain que pour les 
faibles valeurs de r on a une force de répulsion à 


variation rapide, par exemple » $ étant 


(r—r} 
élevé; l’étude exacte de cette répulsion est du 


2: 


y 
e 

- 

" 
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domaine de la mécanique quantique; le cas de molé- 


…cules sphériques impénétrables conduirait à $— 0. 


de 


Il ressort de cette discussion rapide que la loi 


de force est, en général, assez compliquée et mal 
connue; heureusement du point de vue formel il 


SR 5 tn, 


+ 


nous suffira de savoir que la force est toujours cen- 
m…trale; c'est seulement pour préciser les valeurs 
"numériques de certains résultats qu’il est nécessaire 


de connaître exactement la loi de force; les phé- 


… nomènes dépendront d’ailleurs dans une certaine 
mesure de la densité du gaz; en effet, pour les faibles 


Le 


densités les interactions à longue distance sont 


relativement plus importantes que pour les fortes 


densités. 


d. THÉORIES DES GAZ DE LORENTZ ÉLECTRONS- 


- MOLÉCULES EN - PRÉSENCE DE CHAMPS ÉLECTRO- 


” MAGNÉTIQUES. — Si l’on soumet un gaz de Lorentz 


DSNCIPEUS [7 jan 


at cut toi dr ot dd nd à em ve -na dbé sn Colt à Éd RS te db te ile" bin nu à lost on it te + etifé és, CS Le ns De à 


fi duels. … 


- électrons-molécules à un champ électrique E et 
- à un champ magnétique H la distribution des vitesses 


des électrons est différente de celle qu’on aurait en 
l'absence de champs; l’étude théorique de cette 


-distribution et des. propriétés électromagnétiques 


macroscopiques du gaz qui en découlent a été abordée 
par de nombreux auteurs en commençant par les 
cas les plus simples : 

Lorentz [3] a étudié le cas d’un champ électrique Æ 
continu sans champ magnétique; 


_ Chapman et Cowling [4], le cas d’un champ E 
_ continu superposé 


à un champ magnétique H 
continu ; 

Margenau (5) le cas d’un champ Æ alternatif 
sans champ magnétique; 

Jancel et Kahan [6] le cas d’un champ £ alter- 


… natif superposé à un champ magnétique H continu. 


Le cas des champs magnétiques alternatifs est peu 
intéressant car en pratique un champ H alternatif 
est toujours associé d’après les équations de Maxwell 


= à un champ E alternatif et l’ordre de grandeur du 


rapport ?, est tel que les forces magnétiques agissant 


sur les électrons sont négligeables devant les forces 
électriques tant que la vitesse des électrons est faible 
devant la vitesse de la lumière. Le problème traité 


. par Jancel et Kahan représente donc pratiquement 


le cas général. Toutefois, nous avons montré [7] 
que leur formalisme n’est pas adapté à ce cas et 


_ nous y reviendrons plus loin. 


Nous avons été amenés de ce fait à reprendre 


la question de façon plus complète; afin d’établir 


une théorie cohérente et rigoureuse, nous nous 
sommes tout d’abord limité strictement à l’étude 
du gaz de Lorentz parfait; les paragraphes 2, 3, 4, 5 
de cet article et de l’article suivant sont consacrés à 
cette étude; nous montrerons, dans un troisième 


_ article que les hypothèses du gaz de Lorentz par- 


fait ne sont pas exactement réalisées en pratique et 
analyserons les conséquences qui en résultent. 
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2. Exposé de la méthode. — a. ÉQUATION 
DE BOLTZMANN. — Soient f(v, t) la fonction de dis- 


tribution des électrons à l’instant { et F(V, ?) celle 
des molécules neutres; le gaz est supposé uniforme 
et indéfini, c’est-à-dire que f et F ne dépendent pas 
du point considéré et qu’on a à tout instant et en 


tout point | 
J JV) do A, 


fl “, t) dy — N, 


(la notation do représente un élément de volume 


(3) 


(4) 


dv, dv, dv. dans l’espace des vitesses). Les deux 
fonctions f et F satisfont à deux équations de Boltz- 
mann; l'équation de Boltzmann s'écrit de façon 
générale : 

Of of 


Port Bat 


dans laquelle 


“ — gradient de / dans l’espace des positions (ou 
espace de configuration); 
— gradient de f dans l’espace des vitesses 


(ou espace hodographe); 
g, est l’accélération imposée aux particules, c’est- 
à-dire la force par unité de masse. 


Cette équation exprime que la variation au cours 
du temps de f est due aux phénomènes de diffusion 
(—v?) à l’action des forces appliquées te 7) 
et des chocs entre particules représentés par le 
terme C. 

Dans le cas présent, nous avons {pour f et F les 
deux équations : 


0F 


ot 


ot 


=A(F), (5) 


OI J 


+(D+v AR) = J(f). (6) 


Dans les deux équations, nous avons négligé les 


À Q : ; : 
termes en D? le gaz étant supposé uniforme; pour 


D) 
ov 
(pas de force appliquée); le terme de chocs A(F) ne 
représente d’autre part que les chocs entre molé- 
cules;” par suite de la relation (2), les chocs élec- 
tron-molécule sont en effet sans influence sur la 
distribution des vitesses des molécules. En défi- 
nitive, l’équation en F est l’équation habituelle 
des gaz uniformes dont la solution stationnaire 
est la distribution maxwellienne : 


les molécules il n'y a pas non plus de terme en 
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Dans l'équation en f nous avons posé : 


r=°£E, 9 -°4 


: €. 8 
m m Se ho) 


e devant être pris négatif dans le cas des électrons. 
Le terme J(f) représente l’effet des chocs; confor- 
mément aux hypothèses exposées au paragraphe 1 
nous ne tenons compte que des chocs électrons 
molécules. On montre dans les traités de théorie 
cinétique que J(f) est fournie par l’intégrale: 


JC) = jji]] ' F'— fF).gp dp de dy. (9) 


Rappelons le sens de cette intégrale quintuple; 
un choc p—P est caractérisé par la donnée de 
huit paramètres : les composantes D,, v,, 0: et V,, 
V,, V: des deux vecteurs vitesses des particules 
longtemps avant le choc, et deux paramètres géo- 
métriques précisant les positions relatives de ces 
deux vecteurs; les deux paramètres choisis ici 
sont ainsi définis : on considère un plan Il passant 
à tout instant par la molécule P et perpendiculaire 
à la vitesse relative des deux particules avant le 
choc pour un observateur se déplaçant avec ce plan 


E’ 


Fig. 1. 


l’asymptote à la trajectoire de l’électron avant le 
choc est une droite EE’ perpendiculaire en E à I 
(fig. 1); la position de cette asymptote est déter- 
minée par les deux coordonnées polaires p et € du 
point E dans le plan I. 

Ceci dit dans l’intégrale (9) on suppose donné v 
et l’on fait varier V, p et £; on a posé : 


AO 
PINS, 
L—=IV NN |; 


Longtemps après le choc les deux vecteurs vitesses 
sont devenus v’ et V', ces deux vecteurs sont fonc- 
tion de V, pete et l’on a posé dans (9) 


REA 


F'= F(V', t) 


On peut remarquer que de manière générale les 
cinq variables d'intégration figurant dans J ne se 


séparent pas; la forme de J se simplifie, au contraire, 
dans le cas du gaz de Lorentz parfait si l’on peut, 
de plus, négliger V devant v, on a alors : 


VV 


(le choc d’un électron ne perturbe pas le mouvement 
d’une molécule), 


v' est indépendant de V (pour la même raison) 

St (10) 

Dans ces conditions, l'intégration sur V se sépare 
des intégrations sur pete; on a 


NUE J) F(V, t)dv fn vp dp de, 


soit en effectuant la première intégration et en fai- 


sant sortir v 


Er 
ri 


(ff) p dp de; 


72 


HE Ne ff 


J est maintenant un opérateur agissant sur la 
fonction f; cet opérateur est indépendant dei 
est linéaire; par contre, sa forme exacte dépend de la 
loi de choc (à travers f'). 


(11) 


b. DÉVELOPPEMENT DE L'ÉQUATION DE BOLTZ- 
MANN EN FONCTION DÜ PARAMÈTRE L. — Pour 
obtenir des solutions de l’équation (6) on est conduit 
d’abord à développer f en série de fonctions sphé- 
riques; ce choix se verra justifié plus loin. On trouve 
alors que l’équation de Boltzmann est équivalente 
à une suite infinie d'équations différentielles linéaires, 
couplées entre elles par l'intermédiaire du champ 
électrique F. Si l’on suppose que F est faible, on peut 
obtenir des solutions convergentes pour ces équations, 
et nous nous limiterons à ce cas dans cet article. 


L'exposé se simplifie si l’on commence par déve-, 


Nous poserons donc 
F= + fit fat 


lopper f en fonction de F. 
(12) 


x . > { 
en supposant par hypothèse que /, est indépendant 
de Tr, que f, est du premier ordre en F, f; du 
deuxième ordre, 


etc., L n'est pas un paramètre . 


sans dimension; en fait, la suite du calcul montre + 


que pour que le développement de f soit convergent, 
il faut que F soit assez petit pour que l'accroissement 


de vitesse pris par chaque électron entre deux chocs. 


du fait du champ électrique soit faible devant sa 
vitesse moyenne, c’est-à-dire que l’on ait : 
LAN 


— ex, 
Y Up 


* 


s 


(43) 


> étant la fréquence moyenne des chocs que nous À 


définirons ultérieurement. 


Introduisons donc le développement (12) dans. 


Pi 


N<L 
\ 


+ 


équation (9) et séparons les termes d'ordre 6j 1, 274 
Don Lil vient: 


of DR À : 
. 0 M VU) (44) 
+ dfi . mof de se 
| PA Ne. = JG) (45) 
df} dfr dfs NE 
D Re JC) Cs 
E é 
…—. c. DÉVELOPPEMENT DES FONCTIONS DE DISTRI- 


| BUTION EN SÉRIE DE FONCTIONS SPHÉRIQUES. — 
… Pour satisfaire les équations (14), (15), (16) etc., 
nous allons maintenant développer les diverses fonc- 
tions ie ho LE .… en série de fonctions sphériques 
4 ou plus précisément poser : 


Ds” à Er 1 = 
\ fo= D aim in + dyy Sim) (17) 
£ Ji => DCS Cm + bin Sim) (18) 
Ci Lt 14 
È f= D (Gin Cim + Sim), (iIGA) 
ar ur 

F: ; 

n Les fonctions C, et S,, sont les fonctions sphé- 
1 riques classiques 4 
\ Cned O;neosmo, (20) 
F Sin VO mn Sin mo ; (21) 


1 peut varier de o à æ et pour chaque valeur, m peut 
varier de o à { mais Sy — 0; @ est fonction de 0 seu- 
“lement. Le tableau I donne en coordonnées sphé- 
-riques et en coordonnées cartésiennes les expres- 
sions [9] des premières fonctions Cm et Sim. 


_Taszeau [. 


L 
; Coordonnées Coordonnées + 
E sphériques. cartésiennes. 
î PP RERSES ne 1 
3 Cr mieu «RUE v cos Vz 
- CES 3 % Sin D cos ® Dee 
TN ME v sinÜ sin © Uy- 


v? 3 V2 — v? 
Co9 PE fre — (3 cos?) — 1) . 
D » 
Csy....... :3v2sin0 cosb cos o 30-Ve 
ANSE ENT 3v?sin 0 cosÜ sin 30:0y 
MC 3sin20 cos2® 3(v2— vi) 
CRAN RTE 3 sin?0 sin2v Eve vx 


Nous supposons que les coefficients dm», ins. SOnt 
des fonctions « isotropes » de v, c’est-à-dire des 
fonctions de v—|vl; nous devrons d’ailleurs 
supposer qu'ils sont aussi fonction de { pour pouvoir 
- représenter les distributions variables avec le temps 
qui se produisent nécessairement si T est lui-même 
| fonction de {. Ceci dit, si l’on porte les développe- 
. ments (17), (18) et (19) dans les équations (14), (15) 
_et (16) on obtient en effectuant les calculs et ordon- 
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cp 


nant par rapport aux fonetions Sphériques un système 
du type 


Ÿ (aim Cim + cs Sym) = 0, 


(. e L . = 
D (im Cim + Bin SD On (A) 
> (Y4m Cim+ D S Um) ==“) 
qu'on peut satisfaire en imposant 
Lim Un Br B,, VU on 0 ( B) 


On obtient ainsi un système d’équations différen- 
tielles en { déterminant les coefficients @y», Gyms ce. 
Si les coefficients œm, Bim, … étaient indépendants 
de v les systèmes (A) et (B) seraient équivalents à la 
seule condition que les C7, et S;», forment une suite 
de fonctions orthogonales dans un domaine quel- 
conque de l’espace des vitesses. En fait nous avons 
supposé que les coefficients 4», &,, étaient non pas 
des constantes à chaque instant, mais des fonctions 
isotropes de v; cela simplifie beaucoup les calculs; 
on peut s’en rendre compte en remarquant que 4 
représente alors à lui tout seul la fonction de dis- 
tribution isotrope non perturbée. Dans ces condi- 
tions il est facile de voir que les &ym», ay sont aussi 
des fonctions isotropes de v et que les systèmes (A) et 
(B) ne sont alors équivalents que si les fonctions C/, 
et Sy sont othogonales sur une sphère quelconque 
centrée à l'origine de l’espace des vitesses : d’où 
l’emploi des fonctions sphériques. Ces fonctions 


: $ ; D) 
ont d’ailleurs un autre avantage : les opérateurs T REY 
C 


v /\ & . et J qui interviennent dans (14), (15) 


et (16) donnent des résultats relativement simples 
quand ils agissent sur une fonction sphérique comme 
nous allons le voir. 


" he D) 
d. CHOIX DES AXES. ACTION DES OPÉRATEURS T de 
€ 


Ü 
v/\@ es J SUR UNE FONCTION SPHÉRIQUE. — Nous 
pouvons sans restreindre la généralité odu prblème 


supposer que @& est dirigé suivant Oz et que T est 
dans le plan yOz (fig. 2) 


ET 
OP OVER 


1 # 
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Ceci dit, nous avons les règles de calculs suivantes : 


dCi T sind ; 
r TP — = [ STE 
ov 2 
—(l + m)(l+ m1 SE im] 
+ L'eosd(Z+ m) Cham, (29) 
28} T sin 
DE = | CL a,m+1 
ov 2 
+(l+ m)(l+ m1) Crn,mil 
+ Feosb(l + m) Siam; 95. 
OC 
VIN OI mQS 24 
/\ 9 im) ( ) 
l DS Ut We 
VA QE 0 Os (25) 
Ov 
TCin)= = vICins (26) 
JSin)=— VS Im (27 ) 


Ces diverses formules (1) sont démontrées dans 
l’appendice qui paraîtra à la fin de l’article 
suivant; les deux dernières sont tout à fait 
remarquables : elles montrent que les fonctions 
sphériques sont des fonctions propres de l’opé- 
rateur J, la valeur propre v, ne dépendant que de 
l'indice 1. Ceci provient du fait que par suite de 
l'hypothèse des forces centrales J commute avec 
le groupe des rotations; il faut toutefois remarquer 
que dans le calcul de J on considère alors » comme 
un paramètre et que les valeurs propres », sont 
des fonctions de v. On peut aussi remarquer que la 
première d’entre elles », est nulle puisque, d’après 
(10) et (11), J(f) est nul quand f est une fonction 
isotrope. 


3. Calcul de l’approximation d'ordre zéro. — 
a. CALCUL DES COEFFICIENTS Gm. — Portons le 
développement (17) dans l’équation (14); il vient 
en utilisant les formules (24) à (27) : 


Q f + ./ Q ! 
>» C&im Cin + dim Sim) DE dam Cim+ Gin Sn) 


& (2 
= DT (an C'im — jm SP De 
en ordonnant par rapport aux fonctions sphériques 
et en annulant le coefficient de chacune nous obte- 
nons : 


vo = © ( Cow); (28 ) 
io + io = 0 (Cho), (29) 
a Emi 0): (Cia); 

dia QG (St), 

ü, +YIA, =0 (Cr), (30) 
È ! / re 
in + Yläm—= MA Gn (Clim) Lee ae (31) 
5 ; x (MAO) 
im + Wm=— MO Im ( Sim) | ( 32 ) 


(:) Dans l’équation (8) du travail de Jancel et Kahan, 


mn, sein 
le terme en PA est considéré par erreur comme un terme 
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SD NES L'EC L < AE nt LR PACE Es 
$ C \ Ke! 747 RER RES à 
} RU ne Ve 140508 TER 
; os x -"Ÿ ut? 
| VE PNEU IE 


La première de ces équations est singulière; sa 
solution est 
Goo(£) = &o(0)- (33) 
‘Pour les autres, on obtient sans difficulté les 
solutions suivantes : 


al = Aer, (34) 
dim = A {mn TV COS( Eu — mQt), (33) 
Ari = A im et Sin (1m — MO), (36) 


Dans ces formules les quantités A, et ®» sont 
des fonctions de v complètement déterminées si 
l’on se donne la fonction de distribution des vitesses 
à l'instant { — 0. 


b. INTERPRÉTATION PHYSIQUE. — Nous verrons 
au paragraphe suivant que tous les coefficients », 


| 


sont positifs; dans ces conditions on peut faire les 


remarques suivantes 


- Quelle que soit la forme de la distribution f, à 
l'instant { — o elle devient isotrope au bout d’un 
certain temps; les anisotropies initiales repré- 
sentées par les coefficients a» et &, (10) me 
demeurent que pendant une période transitoire. 
La décroissance de chaque ordre d’anisotropies est 
caractérisée par le nombre », qu’on peut appeler 


fréquence de relaxation des anisotropies d'ordre L » 


Divers auteurs avant nous [4], [6], [7] ont considéré 


le nombre », et montré qu'il était égal à la fré- 


quence des chocs dans le cas de molécules impéné- 
trables sans interaction à distance; cette inter- 
prétation est impossible dans le cas général où il 


y a des interactions à distance, car on ne peut plus « 


définir une fréquence de choc. Notre calcul montre 
quelle est la signification profonde des nombres »} 


Quand on a atteint le régime permanent, la fonc- 
tion de distribution isotrope a,,({) qui reste est 
identique à la partie isotrope de la fonction de 
distribution initiale; ceci 
évidente des hypothèses de départ; la formule (10) 
montre en effet que les électrons dont la vitesse 
a une valeur absolue donnée forment une classe 
qui se conserve dans les chocs. 


c. CALCUL DES FRÉQUENCES DE RELAXATION y. 


— Pour expliciter un peu », appliquons la for-, 


mule (26) au cas particulier : 


TOI; VER) 


représenté sur la figure 3. 


en F,v,; c’est en réalité avec nos notations un terme en »,: 
Il s’ensuit une erreur que nous avons signalée [7] dans l’ex- 
pression du tenseur de conductivité ainsi que dans les autres 
formules qui en sont déduites. D’ailleurs la forme que ces 
auteurs ont choisie pour la fonction de distribution n’est pas 
assez générale et ne permet pas de représenter les termes 
en v,. 


est une conséquence | 


La colatitude 0'—7y après le choc n’est autre 
que la déviation; elle ne dépend que de p. 
. Dans ces conditions, on a :. 


Cro(8) = v!Py(cos8) 
2 (P, polynome de Legendre d’ordre Î), 
Co). —=0!, 
Cro(N ) _— vlP(cosy), 
J[Co(0)] = — v:Cw(o) 


2 


= of" ) [out Py(cos y) — v!]p dp de, 


vue on Not [ Tee ro dp, 


0 


san f [1— Py(cos y )] p dp. (37) 
0 


Cette formule généralise celles obtenues par 
. Chapman et Cowling (10) pour 1 = 1 et ! — 2. Si l’on 


À 


Z 


Fig. 3. 


» connaît la loi de force, on peut calculer y en fonction 
- de p et » et en déduire par la formule ci-dessus les 
- diverses fréquences de relaxation. Elles sont tou- 
- jours positives -puisqu'on a toujours P,(cos y) 1. 
Dans le cas des molécules sphériques impénétrables, 
x n ’est fonction que de p, les », sont proportionnels 
à v, ce qui correspond comme on sait à un libre 
parcours moyen indépendant de la vitesse. Dans le 
cas plus général HE loi de forces en r*, v, est 


indépendant de v, résultat remarquable, connu 
depuis Maxwell [11] pour la fréquence », et qui 
simplifie énormément les calculs [12]. 


| 
proportionnel à D PA; dans le cas où s —5, v est 


. 

Ë 

à 

- 4. Calcul de l'approximation du premier 
- ordre. — a. CALCUL DES COEFFICIENTS bym. — Sup- 
. posons maintenant le champ électrique alternatif en 
posant : 


T = y cosut (38) 


ntdit | done GT 
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et portons les développements (17) et (18) dans 


l'équation (15); nous obtenons : 


A . 7 " 
D (bin C'im + bin Sin) 
\'4 oa day 
ne > a nr Sn) 


OCm ’ OSYm 
Dame 
ne (aim nu + mn dv 
Q: / 
na Ÿ V1 ( Dim Cm La TE) 


2 D ( bin Cim— dim Sim )e 


(39) 


Pour simplifier le calcul, nous supposerons l’ins- 
tant { choisi de façon telle que les anisotropies tran- 
sitoires de la fonction f, soient déjà devenues 
négligeables; à partir de cet instant que nous appel- 
lerons {, nous pouvons réécrire (39) sous la forme : 


> (him Cire au bimS 1m) 


Dr: Evi( Dim Cim + Dim Sn) 


Y däw Css GS Tao dCoo 
ov Ov 


= oÿm (Cr — OrnSim)s 


+TŸ 


OT on a 
TvCoo= TC sin bo sin 0 sino + T cosVv cost 
Ù = l'sindS;; + T cost Ci 

et 


dCoo QUE 
ONE 


d’où 
ba ( Dim Cim GE bim Sim) 


4 ADN | 
Gr DUT nl Dim Sim) 


L “ 


à DL sin LS 11 + F cosd C0) 


+ se Cin ER Dim Sim ), 


d’où le système 


boo 0 ‘ (Co); 


1 da 
bo 2 LES ES y cosdcoswé = 0 (C0); 
v dv 
Dan + va bus = Qb1à 2 (Ci 
: 1 da 2 
bla + 1041 + = —Mysinbcoswt——Qb; (Su), 
v dv 
Cv bL = 0 (Ci), 
Dim + VidDIm= MmQ Dim (Cin) 
2/ ! 
dim ST ER Dim = —m Q Dim (Sim) 


(43) 
(4) 


(45) 
(46) 


Il est facile de voir que toutes les anisotropies 
initiales d’ordre ! supérieur à 1 que l’on pourrait 
avoir dans f, s’amortissent comme celles de f,; 
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les équations (44), (45) et (46) conduisent en effet 
à des formules exactement semblables aux for- 
mules (35) et (36); on trouve de même que b,, ne 
varie pas avec le temps; si à l’instant où l’on applique 
le champ le plasma est isotrope et b,, nul on peut 
voir que b,, reste identiquement nul (?). 

D'autre part : 


bio = Bioe at 
( 1 1 da 

R) 00 
lio+viv dv 


bin = Bune-itcos(dis — Qt) 


nf cosdeiw/ b ( A7) 


Q 1 da \ a) 
R jee SEE he QAR rt 
su lGo+w)+Qv dv r:sin- pe ( (48) 
bi = Bne "itsin(bi;— 0?) 
te ne 


a nie l A 
— y sindeiwt 49) 
l{tow+v):+ Ov dv | £ 1 ( 


avec 
R {3} — partie réelle de:z. 


Nous voyons donc que pour les anisotropies 
d'ordre 1 il existe un régime transitoire pendant 
lequel les termes contenant le facteur e-":! s’amor- 
tissent; les amplitudes et les phases de ces termes 
sont déterminées par les conditions initiales. Il 
reste ensuite un régime permanent dans lequel les 
anisotropies du premier ordre sont des fonctions 
sinusoïdales de {; ces fonctions sont proportionnelles 
au champ électrique et ont la même fréquence que 
lui; elles vont nous permettre de définir le tenseur 
de conductivité du milieu. 


b. TENSEUR DE CONDUCTIVITÉ. — Par défi- 
nition, les composantes 7,, 7,, j- du vecteur densité 
de courant j sont les valeurs des courants totaux 
traversant une surface de : em? perpendiculaire 
à l’un des axes, on a donc : 


ae fl Vrf dvx dv, dv;= e Il Cf dvx dv, do:, 
Do ie ff 7 dix dv, dv; = e Î Su f dvxdv, dv, 


Ve ef] v- f dvx dvy dv; = e ff Cu dos do, do. 


Pour faire les intégrations on passe en coordon- 
nées sphériques dans l’espace des vitesses et on 
effectue d’abord les intégrations par rapport à Ü 
et o en laissant v constant; par suite de l’orthogo- 
nalité des fonctions sphériques, tous les termes 
rectangles disparaissent; autrement dit, si nous 


(?) Pour déterminer la valeur de b, lorsque le plasma 
présente des anisotropies initiales, il faut en toute rigueur 
tenir compte exactement de ce qui se passe pendant la période 
transitoire en introduisant dans (39) les expressions complètes 
des coefficients a, données dans les formules (35) et (36) 
et en supposant qu’à l'instant où l’on applique le champ 
électrique la fonction b,, est nulle. 
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appelons do l’élément de surface sur la sphère de l 
centre O, de rayon » nous obtenons : | 


po ve s 
Jette I bas (Cu) dv ds — Le fou vi dv, 


eff Ste die TE four vt do, (50) 


ef] b:0 (C0)? dv ds = Te fo v+ do. 


Pour trouver le tenseur de conductivité sous sa 
forme générale, il nous faut maintenant effectuer 
une rotation d’axes autour de Oz afin d'obtenir 
les composantes de j dans un système d’axes où 
aucune des composantes de T n’est nulle, soient 


Lee 
h 
[l 


Fig. 4. 


OXYZ ce nouveau système (fig. 4). On a : 


LEA a 


= | (ucose — 84 sin ®)v# dv, 


JX=JæCos? — j,Sing = 


+ 


. . : « AT € 
JIY= JxSM9 + Jy;Cos9 — - 


LEE … [ bsovtdo. 


En combinant ces dernières relations avec les 


solutions permanentes tirées de (47) ,(48) et (49) 
et en tenant compte des relations : 


Sin VIT ee VE NON 
Sin Sing —=—Yy,— m Ex 
; e 
+ sinb'coss = Ye— KT, 
: ASE ÿ m 
nous obtenons finalement l'équation complexe : 
ea + 
J=lIsi|é; 
avec 


ë = E, eio!, TUE (3) 


1 Go 
Isl=|—& s 0 (#) 
0 (e) Go 


US Le 


N° 12. 
Fate ee I da : 
Ge) 5 ve, 2 
3m J, (éw + vi) dv 
RENTE AO , d We 
EU. 47 AO) de k de (53) 
3m J, Coeur Oo dv 
re? Ne Q OR L 
D2 = - v° = dv. (54 
“ 3m J (idw+wm}+O2. dv Ge 


Les coefficients diagonaux sont les conductivités 
suivant les trois axes, le coefficient croisé représente 
l'effet Hall. On: voit que : 


— La conductivité 5, dans la direction du champ 
magnétique est indépendante de l'intensité du 
champ magnétique; autrement dit elle est égale 
à la conductivité du milieu en l’absence de champ 
magnétique ; 

=— La conductivité so, dans les deux directions 
perpendiculaires à H est différente de ©, : c’est le 
phénomène de magnétorésistance; pour les faibles 
valeurs de H c’est un effet du second ordre; 

— Le coefficient ©, qui représente l'effet Hall 
est, pour les faibles valeurs de H, proportionnel 
au champ magnétique; c’est un effet du premier 
ordre. Pour aller plus loin dans la discussion des 


formules (52), (53) et (54), transformons-les en 


intégrant par parties : 
oo LV? 


pre? VRANS fe V5 + 
que | sd li | ones) 
3m l LO + y1/0 2 LU) + Vs 


Il est facile de voir que le premier terme est nul : 
pour que le nombre des électrons soit fini, il faut 


T0 = — 


en effet que ‘+ 4Tv’av, dv soit finie, c’est-à-dire 


ue 2&p, décroisse plus vite que r quand - tend vers 
q 00 D 
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l'infini donc que wa,, tende vers zéro; (55) peut 
donc s’écrire 


a = ATv? e?2 I e 
a= f au d 5 (4 
ae TU 
= | an dE cit |. 
à à 
On obtient de même 
SE av d [2 #6) | (87) 
0 
#  ATEUNN, “ 
Go — av d | 3 s, (0) | (58) 


on 


Les quantités ©; (v), s,(v), 2 (v) sont des «conduc- 
tivités par électron » d’expression 


Ë RCE ES : 
FANS F0 HV; ? (59) 
: ee Q 60 
MO Goes et Sue 

2 (eo) 
s(0) = © © (61) 


m (JO +)? + Q? 


Dans le cas où », est indépendant de » on peut 
faire sortir cç(v), si(v), c,(v) des intégrales (56), 
(7) et (58); il reste à calculer 


fa 4Tv°? dv. 


Nous démontrerons dans le prochain article que 
cette intégrale est égale à n, c’est-à-dire que les contri- 
butions à la densité des fonctions f,, fo sont 
nulles. Quand », est constant, on trouve donc : 


Où ai Cie 7001, On = 100, 


Manuscrit reçu le 17 mars 1954. 
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REMARQUES SUR L'ORGANISATION DU CORTÈGE ET DU NOYAU DES ATOMES 
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Sommaire. — On n'avait pu jusqu’à présent introduire la notion de couche dans la théorie des 
systèmes de corpuscules sans négliger une partie importante de la corrélation entre les positions des 
corpuscules. On décrit un procédé qui permet d’associer aux couches d’un système certains domaines 
d'espace. Cette méthode est applicable à la fonction d’onde rigoureuse. Il devient donc possible de 
parler de couche sans négliger aucune partie de 1a dite corrélation. On obtient de plus des ind'cations 
intéressantes sur la localisabilité des corpuscules dans les cortèges et les noyaux qui permettent de 
présenter à partir de modèles géométriques les résultats de la Mécanique ondulatoire sans en heurter 
les principes. Il est intéressant de noter que dans cette théorie plus rigoureuse il semble qu’on sera 
conduit à remplacer les nombres quantiques individuels entiers des théories approximatives par des 
valeurs moyennes seulement approximativement entières, mais nous n’avons pas terminé la discussion 
de cette question qui devrait peut-être ouvrir une voie vers des méthodes de calcul améliorées et 
assouplies des grandeurs observables à partir de la notion de couche. 


Introduction. —— Le succès rencontré, lors de 
l'introduction de la notion de couche dans les 
noyaux d’atomes suggère une analogie intéressante 
entre l’organisation du cortège et celle du noyau 
des atomes. En prévision d’un Colloque réservé 
aux interactions entre cortège et noyau, il nous a 
paru opportun de présenter ici un certain nombre 
d'observations sur la structure respective de ces 
deux systèmes. 

Nous rappellerons d’abord les difficultés que 
l’on rencontre dans l’emploi classique de la notion 
de couche. Plus précisément nous soulignerons la 
nature de l’approximation dans laquelle on est 
conduit à se placer. 

Puis nous montrerons qu’il nous semble exister 
une voie permettant l'introduction de la notion 
de couche dans un cadre plus rigoureux. Nous 
essaierons ainsi d'expliquer pourquoi l'emploi de 
la notion de couche conduit dans le domaine nucléaire 
à des résultats meilleurs que ceux auxquels on 
pourrait s'attendre parce qu’elle entraîne en fait 
moins d’approximations que l’on croit en faire. 
Nous signalerons enfin que l'introduction de la 
notion de couche selon la voie que nous proposons 
entraîne un assouplissement de la théorie qui devrait 
permettre d'obtenir quelques améliorations des 
calculs de Physique nucléaire et une compréhension 
plus claire des phénomènes qui raissent de l’inter- 
. action entre le noyau et son cortège. 


La notion de couche dans le cortège électro- 
nique. — Il est bien connu que d’un point de vue 
expérimental la notion de couche se dégage prin- 
cipalement de l’étude des phénomènes d’ionisation 
des atomes. 

Quand on essaie de retirer successivement les 


électrons d’un atome, on sait que l’énergie minimum 
d’ionisation croît au fur et.à mesure que le nombre 
des électrons restants diminue et surtout que cette 
variation présente des irrégularités (voir, par exemple, 
la figure 1), ce qui donne l'impression que l’on 


5 


a] 
=] 


100 


Energies d'ionisation du bore (en ev) 


ERP 1 
Nombre d'électrons présents 


Fig. 1. 


rencontre dans l'atome des couches d'électrons 
d'autant plus difficiles à extraire que l’on se rapproche 
du noyau. 

On observe des irrégularités dans la variation de 
l'énergie de liaison moyenne et minimum des 
électrons gravitant autour d’un noyau Z donné 
avec le nombre des électrons du système (fig. 2). 

On ne rencontre pas, au contraire, d’irrégularités 
quand on étudie la variation de l’énergie totale de 
liaison des électrons d’un atome neutre dans l’état 
fondamental avec le numéro atomique (fig. 3). 

Il faut se souvenir que l'énergie d’un » niveau 
atomique » est au signe près défini expérimenta- 
lement comme la différence entre l'énergie de 


- liaison d’un atome et l'énergie d’un des ions corres- 


pondants à cet atome : ce n’est donc pas une 
grandeur observable au sens de la Mécanique ondu- 
Jlatoire, mais la différence entre deux grandeurs 
observables. 

Il est encore important de noter que la mesure 
d’une énergie d’ionisation ne fournit pas l'énergie 
d’un électron particulier, mais la différence entre 
l'énergie de liaison de l’ion et celle de l’atome compte 
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tenu de l’énergie de la réorganisation accompagnant 
le départ de l’électron. 
Du point de vue théorique et depuis que l’on a 


abandonné l’idée de trajectoire simple pour les” 


électrons d’un atome on ne sait pas, à partir d’un 
seul état de cet atome, définir rigoureusement des 
niveaux d'énergie pour les électrons. Une seule 
énergie propre a un sens rigoureux : celle du système 
total. 1 

En raison de l’indiscernabilité des électrons ou, 
si l’on préfère, de l’antisymétrie de la fonction 
d'onde, tous les électrons d’un atome ont même 
énergie moyenne, En conséquence, il n’est pas 
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non plus rigoureux d'associer aux électrons des 
nombres quantiques individuels. Seul l’état total 
du système est quantifié. 

La moins mauvaise technique pour introduire 
théoriquement l’idée de couche consiste à remplacer 
la véritable fonction d’onde par une fonction 
approchée ® ayant l’aspect d’un déterminant cons- 
truit sur autant de fonctions œ6»(i) qu'il y a 
d'électrons dans le système, i désignant l’ensemble 
des coordonnées d’espace et de spin associé à l’un 
des électrons. 

On cherche alors les meilleurs ©» et montre 
qu’elles sont solutions de l’équation du champ self 
consistent 


ICPp= Ep?p 


Cette équation ressemble à l’équation de Schrô- 
dinger pour un seul électron mais l’opérateur 4€ 
contient les autres fonctions o et l’analogie n’est 
donc que formelle. Elle suffit cependant pour que 
l’on puisse associer aux différentes fonctions o» 
des nombres quantiques individuels P et une 
énergie €? qui est voisine d’un des « niveaux ato- 
miques » dont la définition expérimentale vient 
d’être rappelée. 

En raison de la forme de la fonction 


o4(1) 81) Pet) 


94(2) 98(2) 
94(3) 98(3) 


le même électron, soit 1, apparaît comme étant sous 
la dépendance simultanée de tous les groupes de 
nombres quantiques approximatifs À, B, C, ... et 
c'est à cet état de fait que l’on associe la notion 
« d'échange » qui suggère bien facilement des idées 
fausses quand on cherche à associer des idées phy- 
siques à cette notion essentiellement formelle et 
mathématique. 

C’est ainsi que dans un cadre approché se fait 
le raccord entre théorie et expérience. Mais en 
examinant la forme de la fonction d’onde approchée 
utilisée D, on voit que ce pont n’a été établi qu’en 
négligeant une partie importante de la corrélation 
existant entre les positions des électrons, c’est-à- 
dire du fait que la probabilité de rencontrer un 
électron en un point donné dépend de la position 
simultanée des autres. 

C’ést surtout la corrélation liée à la répulsion 
coulombienne qui est sacrifiée, la corrélation de 
Pauli étant assez bien représentée par le carac- 
tère antisymétrique du déterminant ®. 

On peut corriger cette approximation en intro- 
duisant l'interaction de configuration, mais la 
notion de couche disparaît alors en même temps; 
ce qui souligne bien son caractère approximatif 
dans les théories actuelles. 


806 


La notion de couche dans le noyau. — L'idée 
de couche se dégage d’un grand nombre de faits 
de la Physique nucléaire et il n’y a pas lieu de les 
rappeler ici. Nous aimerions cependant remarquer 
que c’est l’étude de la variation de l’énergie moyenne 
de liaison des nucléons avec leur nombre (fig. 4) [1] 
qui dans le domaine nucléaire se rapproche le plus 
des études signalées tout à l’heure dans le domaine 
du cortège. On voit, en effet, apparaître des disconti- 
nuités comme dans le cas du cortège et l’on observe 
comme dans ce cas une baisse d'ensemble de cette 
énergie quand s'élève le nombre des corpuscules 
au moins à partir de 62 nucléons. 

Cependant dans le cas nucléaire certaines des 
discontinuités montrent que l’énergie moyenne de 
liaison des nucléons passe par un maximum lorsque 
« la couche correspondante est pleine ». Ce n’est 
pas le cas dans le domaine du cortège. D’autre part, 
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l'énergie moyenne correspondant au noyau le plus stable. 


pour des noyaux de nombre de masse inférieur 
à 62, l’énergie considérée croît avec le nombre de 
corpuscules, ce qui n'arrive jamais dans le cas du 
cortège. 

Ces différences entre les deux systèmes de cor- 
puscules résultent du fait qu'entre deux nucléons 
l'interaction contient une partie attractive alors 
qu'entre électrons corticaux, elle est toujours 
répulsive et que dans le noyau il n’y a pas de centre 
de force sensiblement fixe comme dans le cas du 
cortège. La stabilité énergétique spéciale aux 
couches nucléaires conduit de plus à penser que 
l'interaction attractive entre les nucléons d’une 
même couche pleine est particulièrement forte. 

Du point de vue théorique, les approximations 
que l’on introduit habituellement pour retrouver 
l’idée de couche sont ici de même nature que dans 
le domaine du cortège. On introduit la corrélation de 
Pauli, l'esprit du champ self et on néglige assez lar- 
gement la corrélation liée au champ spécifiquement 
nucléaire. L'emploi des nombres quantiques indi- 
viduels et la notion « d'échange » présentent encore 
le même caractère équivoque que dans le cas du 
cortège. 

Cependant comme on admet habituellement que 
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la corrélation liée au champ nucléaire est plus impor- 
tante entre nucléons que la corrélation coulom- 
bienne entre électrons, on est en droit de s'attendre 
à ce que les modèles en couches conduisent à des 
résultats assez décevants, ce qui n’est pas tellement 
le cas. 


Théorie de la localisabilité des électrons 
dans le cortège électronique [2]. — a. LA 
NOTION DE LOGE D'ORDRE 1. — HEnvisageons un 
système d’électrons, par exemple. Supposons qu'il 
soit dans un état à projection de spin bien définie 
et soit p le nombre des électrons de spin += Sa: 
par exemple. Découpons l’espace en p volumes V; 
connexes ne se recouvrant pas et recouvrant tout 
l’espace et tels que l’on puisse évaluer dans chacun 


d’eux la probabilité P de rencontrer un électron 
h 


et un seul de spin += 


2 27% 


Nous dirons que l’on a découpé l’espace en loges 
; S JA) 220000 
d'ordre 1 et de caractère + - -— vis-à-vis de.ce 


système. La même définition se transpose au 


ê 1 À 
Spin — Fe SA 


La quantité 
9 D P; 


est supérieurement bornée par la valeur p. Posons 


ses PE Le 
P 
n est positif et nous lui donnerons le nom de défaut 
de localisation du découpage par rapport au système 
dans l’état envisagé. 

Le découpage sera considéré d’autant meilleur 
que n est petit. 

L'ensemble des découpages en loges ainsi définies 
est très vaste, mais dans la pratique nous nous 
bornerons souvent à considérer des sous-ensembles 
en imposant aux découpages des propriétés supplé- 
mentaires suggérées, par exemple, par la symétrie 
de l’état étudié. Il deviendra alors généralement 
possible de trouver le (ou les) meilleur (s) décou- 
page (s) d’un tel sous-ensemble, c’est-à-dire celui 
(ou ceux) correspondant au minimum de quand 
il existe. 

Soit Wir (LG m0 OS 2005 LOS RO TERRE 
la fonction d’onde représentant le système dans 
l’état considéré, on a 


Pi= CP du f ds. tp f. dept ee dEn 
LE 


“ espace-7'; espace 
dis 


" 
do» | R do p+1 ... du, W*w, 
. AVE 


Er 227 


Il est donc facile de calculer les P, correspondant 
à un jeu donné de V;, donc d'évaluer n. 


pute REMARQUES SUR L'ORGANISATION DU 


b. RECHERCHE DES MEILLEURS DÉCOUPAGES EN 
LOGE DANS LE CAS DE L'HÉLIUM. — Discutons de 
ce point de vue deux états excités de l’hélium. 
Envisageons d’abord le premier état triplet. 

Il est bien représenté par la fonction [3] 


Br) —SerZ's short, 
avec 
Ti Fo 


SE —— » 
(220) (421) 


AN 374, c — 0,895, 


Adoptons le cas où les deux électrons ont le 
D, 1e 7 ue s MUR 
IS et La symétrie du problème suggère ici 
deux découpages distinets : soit un plan passant 
par le noyau, soit une sphère de rayon R concen- 
:trique au noyau. Dans les deux cas, on définit deux 
loges où 
Pa Rs tre 


Dans le cas du découpage plan on trouve 


Pi 0,5; 1000? 
P 
1,00 
(a) 
0,75 
0,50 
025 
NAN RES) 4 5 


R (en unités atomiques) 


Fig. 5. 


La courbe a de la figure 5 montre la variation de P 
avec R pour le découpage sphérique. Il existe pour P 
un maximum (P =o,93 pour R-1,7 45), donc 
pour » un minimum (0,07). 

Ce minimum étant plus petit que dans le cas du 
découpage plan, il s'ensuit que de tous les décou- 
pages considérés, la sphère de rayon R —1,74, 
donne le meilleur. 

Notons bien la signification physique de ce 
résultat : il y a 93 chances sur 100 que sur son premier 
état triplet l'hélium comprenne un électron dans la 
dite sphère et l'autre à l’extérieur. 

Envisageons maintenant l’état triplet de pro- 


jection de spin + 6 représenté par la configu- 


ration antisymétrisée (2s)(2p-). Adoptons des 
fonctions hydrogénoïdes avec Z — 1,65. On trouve 
pour le découpage plan (perpendiculaire à l'axe 
des 2) 


P=0,78 et T== 1090: 


: . 
Euh 
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La courbe b de la figure 5 montre que dans le cas 
du découpage radial la grandeur P ne dépasse 
pas 0,52. [ci c’est donc le découpage plan qui est 
le meilleur. 

Dans le cas du premier état triplet qui correspond 
classiquement à un électron X et à un électron L, 
nous avons été conduits à partager l’espace à l’aide 
d’une sphère et l’on peut associer l’espace intérieur 
à cette sphère à la couche X et l’espace extérieur, 
à la couche L. Dans le second état correspondant 
à deux électrons L nous avons été conduits à couper 
l’espace par un plan. Ce fait est en accord avec 
l'existence d’une seule couche qui se sépare en 
deux loges. Dans le premier cas, la corrélation 
radiale prédomine. Dans le second, c’est la corré- 
lation angulaire. Il y a là un moyen de définir 
la notion de couche à partir de n’importe quelle 
fonction d’onde et notamment à partir de la fonc- 
tion d’onde exacte sans faire, bien entendu, usage 
de nombres quantiques individuels. 


€. COURONNES  SPHÉRIQUES  ASSOCIÉES AUX 
COUCHES. — En utilisant ce procédé, on est finalement 
conduit à faire correspondre à la notion de couche 
d'un atome un ensemble de loges qui pavent l’espace 
d'une couronne sphérique concentrique au noyau. 
Les sphères successives qui définissent ces couronnes 
ont [4] sensiblement comme rayons ceux qui corres- 
pondent aux minima successifs de la densité radiale : 


Da)=T | Er AD ONE, 5, 2) 8r ri) dt 
ni 


dans le cas des états à symétrie sphérique. 

La charge électrique moyenne contenue dans la 
couronne ainsi associée à une couche qui aurait 
classiquement le nombre quantique n est voisin 
de 2 n? pour un état à « couches complètes ». 

On a pu (!) sur quelques cas vérifier que l'énergie 
moyenne d’un électron dans la couronne ainsi 
associée à une couche est voisine de la valeur de ep 
que la théorie du champ self lui ferait correspondre. 

On a donc l’impression qu'aux nombres quantiques 
individuels et entiers qu'une théorie approximative 
fait correspondre à la notion de couche il y a lieu de 
faire correspondre des valeurs moyennes approxima- 
livement entières dans une théorie rigoureuse. 

C’est là où se trouve à notre avis la source d’un 
assouplissement qu’il serait souhaitable d'introduire 
dans la notion de couche. Nous ne prétendons pas 
avoir trouvé la voie la meilleure pour l’obtenir, 
nous croyons seulement avoir indiqué une direction 
possible de recherche. 

L'amélioration essentielle que nous obtenons en 
essayant de relier la notion de couche à certains 
domaines d’espaces, à certaines valeurs moyennes 
c'est la possibilité d'utiliser une fonction d’onde 


(2) Calculs inédits. 
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exacte (ce qui n’est pas le cas de la voie classique que l’on peut adopter pour caractériser les diffé- 
comme nous l'avons rappelé). Il n’est donc plus rentes couronnes sphériques associées aux couches, 
nécessaire de négliger une partie de la corrélation dans l’état fondamental de quelques atomes ou 
entre les corpuscules et nous voyons là un avantage ions, d’après les fonctions d'onde du champ seli=M…} 
appréciable notamment dans le cas du noyau. consistent, ainsi que la charge moyenne que l’on 
Le tableau ci-dessous [5] donne les rayons r, rencontre dans chacune d'elles. 


Éléments... Be. Re AIT. Care: Rb+. Hg. 
RÉ EURE 1,125 + 0,00 0,379 + 0,009 0,22 + 0,009 0,13 +0,01 0,06 —+0,0025 0,025 + 0,000 
Charge....,. 9,05 +0,01 2,240 0702 2,21 + 0,04 D, 2 10) 231 HO 2 + 0,05 
ACER - - _- 0,64 +0,01 0,265 + 0,005 0,102 + 0,001 
Charsete - = - SE OS ns 10,29 9,9 MER 
RU Le ES - - - : 1,095 + 0,01 0,278 + 0,005 
Charge... ... - 2 - = 18,10 0,06 18,10 
NES ee - - - a = 0,93 0,09 
Charge... ... - = - - _ 33 HU 


Les valeurs 7,, sont données en unités atomiques, les charges en électrons. 


Ce tableau permet de calculer sans peine le volume Il existe ainsi une loi évoquant un peu la loi de 
associé aux différentes couches et celui V des loges Mariotte entre le volume moyen occupé par un électron 
associées aux électrons dans ces couches. Le volume V ef le potentiel électrique qui agit sur lui. 
donne une idée de « l’espace vital moyen » réservé 
à chaque électron quand il passe dans le domaine 


HAS La notion de loge dans les noyaux. — Nous 
associé à une couche. 


RE à ONE F ___ avons rappelé que la plus grosse objection que l’on 
Si l’on évalue le potentiel électrique moyen P peut présenter contre les modèlés nucléaires « en 
couches » est le fait qu’ils semblent conduire à une 
négligence d’une grande partie des corrélations 
d’origine spécifiquement nucléaire. $ 
Dans le cadre du langage des loges, il devient 
possible de parler de couche tout en conservant la 
totalité de la corrélation puisque l’on pourrait en 
principe définir des loges à partir de la fonction 
d'onde exacte caractérisant un noyau. 


En fait comme on est encore loin d’avoir une 
bonne expression de cette fonction d’onde et comme “ 
une partie de nos efforts a été dirigée vers l’appli-« 
cation prometteuse de l’idée de loges à l’étude de 
la structure des molécules, nous n’avons pas à cette 
date effectué de calculs permettant d’avoir une idée 
un peu précise de l’organisation des loges dans les 
noyaux. Nous nous bornerons donc à quelques 
considérations générales. 


Remarquons d’abord que l’on admet le plus 


& Couche K souvent que le potentiel d'interaction nucléaire est 
sa $ presque constant à l’intérieur du noyau puis diminue 
AR A rapidement en valeur absolue quand on s'approche « 


de la périphérie. Si l’on admet encore ici une loi en: 


ms P°V = const., 


log me Es 
Fig. 6. la constance de P entraîne l'existence de loges 
ayant un volume sensiblement constant, sauf près … 


ressenti par l’électron considéré dans ledit domaine, de la périphérie où ce volume doit devenir plus 


on s'aperçoit que l’on a sensiblement [6] grand. 
L Cette quasi constance du volume moyen occupé … 
Pre ainsi par un nucléon est en bon accord avec le fait 
que le volume du noyau est sensiblement propor- 
(comme le souligne la figure 6). tionnel au nombre des nucléons qu’il contient. 


, 


; 
] 
È 


ZOPE EX, VU NS OS), SET 


RAT À IE TS 4 Le € 


“types de loges : 


19 Les loges pour protons de spin + =; 


20 Les loges pour protons de spin — : ; 


39 Les loges pour neutrons de spin + = 


4° Les loges pour neutrons de spin — 5 

Les loges de types différents peuvent se super- 
poser et l’on peut donc concevoir des régions de 
l’espace où il y a simultanément place pour deux 
protons et deux neutrons de spin contraire suggé- 
rant l’idée de particules « «préformées ». Mais notons 
bien que cette notion se retrouve ici sous une forme 
très assouplie par rapport à celle que l’on rencontre 
dans le cadre de la théorie du modèle à particules & 
préformées. 

Ici l’organisation en particules « n'apparaît que 
d’une façon éphémère d’autant plus que le défaut 
de localisation des meilleurs découpages en loges 
nucléaires est certainement bien plus grand que 
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Remarquons qu’on devra ici considérer quatre 
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dans le cas du cortège électronique (?), c’est-à-dire 
que les loges ne correspondent ici qu’à une locali- 
sation bien moins forte que dans le cas du cortège. 

Finalement, notons encore que la répulsion cou- 
lombienne doit rapidement tendre à augmenter le 
volume des loges pour protons par rapport à celui 
des loges pour neutrons. On comprend donc d’une 
autre façon pourquoi dans un même espace nucléaire 
il y a généralement place pour plus de neutrons 
que de protons. 


Conclusion. — On sent finalement que la notion 
de loge apporte à la fois de la souplesse dans la 


notion de couche et la possibilité de « visualiser » 


sans manque de rigueur les résultats de la Méca- 
nique ondulatoire. Cette possibilité de faire reposer à 
nouveau l'esprit sur des modèles géométriques ne 
heurtant pas les conceptions de la Mécanique ondu- 
latoire nous semble l’un des résultats les plus sédui- 
sants de cette théorie. 


(2) On se rend facilement compte de cet état de chose en 
réfléchissant sur les fonctions de l’oscillateur harmonique 
qui sont souvent utilisées dans la théorie des modèles en 
couches. 


Manuscrit reçu le 10 avril 1954. 
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I. Introduction. 


Nous réunissons 1ci, et nous systématisons, un cer- 
tain nombre de Notes ou de Mémoires antérieurs. 
Il ressort des formules même que nous produisons que 
la symétrie et la covariance relativistes peuvent être 
introduites en Théorie des Quanta dès le niveau élé- 
mentaire de la théorie de la particule à spin non 
explicitement superquantifiée. Dans une autre partie 
de notre étude, publiée ailleurs [9], nous montrons 
que la prise en considération des statistiques quan- 
tiques est indispensable à la levée de certains para- 
doxes, en sorte que, finalement, la superquantifi- 
cation est aussi indispensable que le spin à la 
réconciliation complète de la Relativité et des 
Quanta. Quoi qu’il en soit, nos formules très simples, 
ici produites, viennent corroborer les beaux travaux 
de Tomonaga [18], Schwinger [16], [17], Dyson [11], 
Feynman [12], auxquels elles constituent en somme 
une introduction élémentaire. 

Mais cette étude technique en appelle une autre, 
d'ordre épistémologique, dont le titre tout indiqué 
sera Complémentarité et Relativité [10]. Nous y 
discutons, notamment, en termes d’Univers de 
Minkowski, la célèbre expérience de pensée d’Eins- 
tein, Podolsky, Rosen, et quelques autres. 

Finalement, nous croyons pouvoir affirmer que 
la réconciliation des formalismes de la Relativité 
restreinte, et de la Théorie des Quanta fondée sur la 
Complémentarité, est un fait absolument accompli. 


II. — Particule libre à spin. 


I1.1. Transformations de Fourier cova- 
riantes [7]. — Soient À et c les constantes univer- 
selles bien connues, k; la /{-fréquence d’une onde 
plane et p, l'impulsion-énergie du corpuscule 
associé, k, et cm, les longueurs de ces /4-vecteurs : 


hkh=27p), 


hko—= 27cm, l 


(1) 
TA DNS AD \ 


L’équation de la particule de spin multiple de en 
est de la forme 


adi+ ko} d(x) = 0, (2) 


et, quelle que soit l’algèbre des matrices a, l’équa- 
tion de Gordon 

{di — kB} V(x)= 0 (3) 
est conséquence de la théorie. 


Avec Marcel Riesz [15], prenons le développement 
de Fourier du % sous la forme ( 


kr? (x) = 


appliquant l'opérateur de Gorden, 
CI si AZ KE, 
H=— 2 si (4) Zo: 


et E(k) dix: (4) 
on voit que 
Ainsi, # désignant l’hyperboloïde à deux nappes 


du {-espace k 
Fi + k2 = 0 (5) 


et 


idu=[dAsdkdk,],  Hôv=— Hô, (6) 


définissant le {-vecteur élément de volume sur 4 
(évidemment colinéaire à k;,) ainsi que son module 
(positif par définition), le développement de Fourier 
le plus général du 4 peut être donné sous les deux 
formes équivalentes 


a 
3 


(27) V(æ)=— à ff ser etse (x) 


de 


= ff] Sveret(x): (7) 


àe 


on passe d’une forme à l’autre au moyen de (6;) 


et (2). Naturellement, la forme en k de (2) est 
{ a A — 10 } É(Æ) = 0: (8) 


Je dis que la formule de Fourier covariante réci- 
proque de (7) est 


( 27) t(E) = — JT sue kje—ikha, ( ku + dH) V(æ), (9) 
20 & 

où & désigne une hypersurface quelconque du genre 
espace de l’espace-temps (?), du, son 4-vecteur 


(:) Le début de cette étude s'inspire de Màrcel Riesz [15], 
que nous généralisons sur les points qui nous intéressent. 
(?) Nous préférons la notation & au + de Tomonaga- 


: 
{ 


élément de volume (idu, > o), € 


à 
7 


A. 


jo 42, - 


_ de signe bien connu : 


1 


ki 
A 


FT 


Prenant l'intégrale (9) sur deux & 


tou, = [dry dx, dr, |, E(Æ) — 


(10) 


différentes, 


“transformant la différence en intégrale quadruple, 


et tenant compte de l'équation de Gordon sous ses 


“formes (3) et (5), il vient bien zéro (*) en vertu de 


| MEAIDIErS 


| 


dpt (iku+ QU) (x) | 
= exthe Ki kb— ki du ik ou+ où } 


pa). 
Prenant alors pour & l’hyperplan x, —o, et 
écrivant séparément le résultat pour les énergies 
| positives et négatives, avec 
ne. 


1, 2, 3), (11) 


il vient. 


7 Mae Rë p, + d:)4(%, 0). (42) 
0 


: Par ailleurs, faisant dans (6,) À —4, il apparaît 


linvariant relativiste bien connu de la théorie des 
champs 


3 
Le LT + sur d€,, — sur ÿ€_, (13) 
k 


-de sorte que l’on déduit de (7,) les expressions 


suivantes pour les contributions à - 4 des 


énergies io et Fes 


ou, équivalemment, 


na ue (&, 0) = ko Jak eñ 3e.) (15) 


D'après la théorie élémentaire des intégrales de 


Fourier. les réciproques des (14) et (15) sont 


Gae()=+ E ffañe &rine( 0) 
to 
= + Fe J d% ef 0.(%, o), (16) 
0 


et il s’ensuit bien (12) ci-dessus, si l’on remarque 
que toute composante plane à énergie positive 
du % est orthogonale (au sens ordinaire) à toute 


composante à énergie négative; il y aurait doute 


D “ 
pour deux composantes de même X, mais, comme 
cette égalité n’est que relative, on lève le doute 
en changeant de repère galiléen. 


Schwinger, notamment par ce qu’elle laisse disponible l’habi- 
tuelle notation de la densité de spin. 

(#) Ce calcul, et la formule (9), rectifient ceux que nous 
avions initialement donnés [4], [6]. 
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(k) un commutateur 
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Dans (9) figure la dérivée normale sur & du 4 (x) : 
celà tient à ce que nous n’avons invoqué que l’équa- 
tion du second ordre de Gordon. En faisant jouer 
léquation du premier ordre de la particule à spin, 
on fera disparaître la dérivée normale, mais, en 
contre-partie, l'algèbre particulière des a; fera 
son apparition. Par exemple, en théorie de Dirac 


du Ÿ = | LYvy.] dY— Ko Yu (UE | 
Co Ho 


Mat ( 0 Siu—Y | 
en théorie de Kemmer [13, éq. (6) et (62)]. 
D Y = ET BV Bu — Bu Bv1d)— Xo Bu } Ÿ. (18) 


Portant dans (9) ces expressions de 04, intégrant 
par parties, et remarquant qu'en vertu de l’anti- 
symétrie de la matrice [ ] ci-dessus, la partie 
« tout-intégrée » est effectivement la transformée 
d’une intégrale double prise sur le contour à l'infini 
de &, intégrale qui est nulle sous les hypothèses 
habituelles concernant le comportement du d à 
l'infini spatial, il reste une intégrale où le dépla- 
cement de l’opérateur d” sur l’exponentielle donne 
de nouveaux termes en X'; réductions faites, il 
vient, en théorie de Dirac, 


(A7) 


et, 


“af ,( 


re [eue e—#khæA | CEA Yu ku— ko) y j Q L(x) (19) 
TEA 
et, en théorie de Kemmer, 
(22 Ÿ € (A) BU DIRTORE 
; 24 & ; 
X {RY + TR [BU ÉY — PV EU] — A7} Lx). (20) 


Le lecteur pourra vérifier à nouveau qu’en prenant 
les intégrales (19) ou (20) sur deux & différentes, 
transformant la différence en intégrale quadruple. 
et faisant jouer cette fois l'équation du premier 
ordre de la particule à spin, il vient bien zéro : 
le €(k) défini par (9), (19) ou (20) est indépendant 
de &. 

Dans (19) et (20), la dérivée normale du 4 a 
disparu. À la place, il y a une combinaison linéaire 
des composantes du 4. 


11.2. Égalité de Parseval covariante, — Intro- 


duisant les adjoints covariants % et ? du L et du £, 
ainsi que la matrice 8 que sait définir toute théorie 
de particule à spin, telle que 


—: + 6, 
les adjointes relativistes de (2) et de (8) sont 
OCK) Lo AX—GTk0 


VE — L+p p, (21 ) 


(x) { ay — ko l = 0, | — 0; 
Lg ) 


(8) et (22 


GG = ko ax 6, 


+) ont les conséquences 


Ne GEL RO; 2 2 
Ah a = 1ko616s, 


(23) 
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qui montrent notamment que les 4-vecteurs de 
l’espace des k, k, et Ga@,, sont colinéaires. 
Je dis que la formule de Parseval covariante 


este) 
à ff] dc vaou =1 ff] (RG at:50, 
& àe 
= /} E(k)Ca Gr dv; (24) 
de 


on passe de la première à la seconde forme au moyen 
de (23,) et (6.). Compte tenu de l'équation de 
continuité 

dy Di Ve = 0, (25) 


l’on voit que le premier membre est indépendant 
de &, ce qui nous permet de le calculer en particuliér 
sur l’hyperplan x, — 0. Par aïlleurs, exprimant le 
troisième membre en tenant compte de (13), puis 
de (23,) pour À — 4, il vient (avec le même signe 
pour les deux nappes) 


x 2e 1 F2 Æ S- 
J] Drag: d%=— & f]] 10 
PES à ; ky—0 ki 

\ 


0 


séparant, en vertu d’une précédente remarque, les 
parties à énergies positives et négatives du Ÿ, l’on 
peut écrire, le signe étant partout + ou —, 


x = A PUUE  M ce 
F ; a s 
I Plata drz=— #4 ff LE  - PE GA (27) 
Am 0 


Led mt) 14 


ceci est bien l’égalité de Parseval ordinaire associée 
aux transformations de Fourier ordinaires (14) 
et (16). CE 


11.3. Définition covariante de l’orthogonalité 
et de la normalisation [2]. — Deux 4 ou deux £ 
seront dits orthogonaux au nouveau sens covariant 
si l'intégrale (24) vaut zéro. Prenant 


Lu 
PLEINS 


(24) fournit la norme de chaque fonction d’un sys- 
tème orthonormal (normalisation en nombre de 
particules). Il suit de ce qui précède que l’ensemble 
des composantes planes à énergies (ou fréquences) 
positives et négatives forme un système orthogonal 


complet, au nouveau sens covariant, pour déve- . 


lopper le . 

En théorie de Dirac, dy#4 = d+4; il y a donc 
identité entre la nouvelle définition et la définition 
classique. Remarquons aussi qu’en théorie de 
Dirac le premier membre de (24) est défini positif, 
en sorte qu’au troisième membre €£ est positif 
sur la nappe ä., négatif sur la nappe 4€; nous 
retrouvons ainsi un résultat connu. 


(*) Cette formule rectifie, par le « (k), celle que nous avions 
initialement donnée [4], [6]. 


Cri 


En théorie des particules à spin supérieur à 274 
lintégrale (26,) s'écrit à 


Î] vieavars: 


ceci est bien connu, en ce qui concerne la norma- 
lisation tout au moins (normalisation en nombre 
de particules). Le premier membre de (24) n’est 
plus alors défini positif, en sorte qu’on ne peut 


rien dire touchant le signe de &£ sur chaque nappe. 
de äC. Par contre, il reste vrai que, pour la norma-" 
lisation en nombre de particules, l’on ne change rien 
en permutant simultanément le signe du {-vecteur k:; 
d’une composante spectrale, et celui de son nombre 
d'occupation mx. 

Vérifions directement, sur la première intégrale (24), 
qu’en permutant les indices 1 et 2, on passe au nombre 
complexe conjugué. Les conjuguées des quatre 
composantes de du sont +du,, + si À =1, 2, 3, 
— si À —4. En supposant, ce qui est toujours 
possible, les a* et fB hermitiennes, les conjuguées 
des quatre composantes de Vi Bad, sont 


DEP = dibalds, — si À =1, 2,3, PS4 


(en théorie de Kemmer, voir [13], éq. (3). Fina- 


lement, | 
(ad ah du) )* = dde ab du. (28) 
c, 


Ici, comme en d’autres questions similaires, c’est 
donc essentiellement le jeu de la « matrice-lampon » 8 
qui concilie les exigences des formalismes relativiste et 
quantique, et notamment les rôles des deux symboles i 
propres à chaque théorie. Signalons que ce jeu de 
formules (ou d’autres très similaires) sont impliqués 
à plusieurs reprises dans Schwinger [16, éq. (1.57) 
et 1(2229)|: 


I1.4. Fonction caractéristique de l’impulsion- 
énergie. — Il est clair que, dans (24), l’intégrand. 
des seconds membres est la fonction de distribution 
de l’impulsion-énergie. Posons alors 


h=t L=etnt; (29) 
{ 


compte tenu de (7), (24) devient 
ie# "brad T+y)ou = 1 A]. e(kjei Matt, (30 
URI Jar o( y )0u RG : $ > (30) 


et ceci montre, en vertu d’une définition bien connue 
du Calcul des Probabilités, que le premier membre 
de (30) est la fonction caractéristique de l’impulsion= 
énergie de la particule libre. Nous avions antérieu= 
rement déduit ce résultat autrement [5]. 


11.5. Résolution du problème de Cauchy. — 
Portant (9), (19) ou (20) dans (7), on résout le 
problème de Cauchy sous la forme 


ya) = f] S\(æ — 2')U(x') du, (æ'}, (31) 
… 6". 


{ 
e 


‘avec respectivement, pour la particule de spin 
quelconque, et D (x) désignant la fonction de propa- 
gation de Stueckelberg (réelle; impaire, et nulle 
- dans l’ailleurs) 


LÀ Si(x) = dd D(x) + D(æ)0}, 

| Le Ra M seras, (32) 
ES 

k pour la particule de Dirac 

$ S(1)(æ) = {Cradë— ko)7}} D(&) (33) 
Let pour la particule de Kemmer 

_. Si(r)= (+, [BE PEUT — of} D(x). (34) 


ê, Avec la forme générale (32), impliquant seule- 
ment l'équation du second ordre de Gordon, la 
dérivée normale sur &’ du Ÿ figure. Avec les formes 
poccansées (33) ou (34), l’on a à la place une 
combinaison linéaire des composantes du 4%. (32) 
«a été donnée par Schwinger pour le photon, et (33) 
pour l’électron [16, éq. (2.22), (2.23), (2.24)], 
MII7, éq. (À.29)|. 

i 


11.6. Localisations patio-temporelle des parti- 
cules. — Le problème de la distribution des parti- 
cules sur les états d’impulsion-énergie a été résolu 
aux paragraphes II.1 et 11.2. L'introduction des 
“fonctions de propagation D(x—x') ou S\(x—x') 
Lva nous permettre de traiter le problème complé- 
mentaire de leur localisation spatio-temporelle. 


n décrive une hypersurface arbitraire € du genre 
“espace, ne coupant pas &’, et que & et 6’ soient 
…_pavées d’une infinité de petites cellules jointives; 
“si l’on aime un langage imagé, disons que & et €’ 
«représentent, dans l'Univers de Minkowski, deux 
“écrans tridimensionnels, chacun percé d’une petite 
‘ouverture (l’une des précédentes cellules) que 
“franchit ou ne franchit pas la particule. En termes 
“de théorie de la mesure quantique de Von Neumann, 
la particule répond ainsi « oui» ou «non» à la 
mesure & ou &/ assimilée à une question; idéalement, 
lon dispose de deux jeux complets &' et & d'écrans 
complémentaires (au sens de l'optique), ce qui 
permet de poser de toutes les manières possibles 
“la double question « &' puis & », ou « & puis &' ». 

Par (31), L(x), considérée comme une variété 
attachée à &, est développée sur le système des 
fonctions SÀ (x — x’), où x est considérée comme la 
variable et x’, attachée à &' “6, comme l'indice 
» de numération. L'on voit sans peine que deux S* (x) 
différentes, c’est-à-dire attachées à deux x’ diffé- 
-rents, sont orthogonales sur & au sens covariant 
- du paragraphe 11.3; il suffit pour celà d’exprimer 
l'intégrale curviligne en termes de sa projection à 
- temps constant, et de noter la présence d’un facteur 
- exponentiel oscillant à exposant non nul. Les 

coefficients du développement du % (x) lié à & sont 


. 
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Imaginons que l’instant-point x de la formule (31) 
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les d (x') liés à &’. Reste à voir si nous pouvons 
calculer ces coefficients par la formule habituelle 
dûment généralisée. 

Commençons par le cas de l’électron de Dirac, 
où le formalisme est d’un maniement particuliè- 
rement aisé. Modifiant légèrement les notations, 
nous écrirons, avec Schwinger [16] 


Sa y(2)= Crudt— he) Da) (35) 


2) 
et 


Ya) = fe SCœ— a) rp(a) Gta’), (86) 
ce qui est une expression du type canonique (24). 
Réciproquement, nous aurons la même formule 
avec permutation de x et x’, et substitution de €’ 
à &, ou encore, en prenant l'expression conjuguée 
et multipliant à droite par y{ 


Ja) = ff rate) Ya) SCe — 2); (37) 


ceci est exactement la formule classique pour le 
calcul des coefficients (x) du développement 
du d(x) sur les fonctions S (x — x), compte tenu 
de la généralisation des définitions impliquée dans 
la formule (24). Nous disposons donc d’un formalisme 
covariant relativiste, transposant textuellement le 
formalisme quantique habituel, qui nous permet 
de traiter le problème des mesures de localisation 
successives d’une particule, assimilées au franchis- 
sement d’un écran tridimensionnel du genre espace 
barrant tout l’espace-temps, et percé d’une petite 
ouverture arbitrairement découpée. Dans le fran- 
chissement ou le non-franchissement des écrans 
successifs (réponse oui ou non à la question posée), 
l’on assiste à la réduction du paquet de probabilités, 
bien connue. La différence essentiellé avec le forma- 
lisme « à temps constant » de la classique Théorie 
des Quanta, est que les fonctions de propagation 


de Stueckelberg remplacent les (x) : si une particule 
a été localisée dans une petite ouverture au fran- 
chissement d’un écran & (localisation tridimen- 
sionnelle, ou localisation de franchissement) elle 
est, dans l’Univers, attachée à la fonction de propa- 
gation D (x) ayant son sommet au point moyen 
de l’ouverture, c’est-à-dire que toutes ses manifes- 
tations passées et futures sont contenues à l’in- 
térieur du cône isotrope de sommet x (localisation 
{-dimensionnelle). 

Avec le formalisme de Kemmer, les formules 
seraient un peu plus compliquées. En particulier, 
les BG} n'ayant pas d’inverses, on ne verrait pas 
un $? prendre la place du y} dans les formules 
analogues à (36) et (37); par voie de conséquence, 
ceci impliquerait une légère extension aux défi- 
nitions posées au paragraphe II.3, extension que 
nous laisserons ici de côté. 


II1.7. Cas d’un écran du genre temps. — 
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Nous avons déjà discuté cet intéressant pro- 
blème [3], [7], et nous le faisons de manière plus 
complète dans un article à paraître [9]. Cette expé- 
rience de pensée constitue un moyen très parlant 
d'introduire l'interprétation des énergies négatives 
de Feynman. De plus, sa discussion montre que 
l'intervention de l’une ou l’autre des deux statis- 
tiques quantiques est indispensable pour écarter le 
paradoxe d’une causalité rétrograde, qui appa- 
raitraît sans celà; le paradoxe en question se trouve 
écarté par le trait, nouveau, des statistiques quan- 
tiques, qui.est de faire intervenir symétriquement 
les nombres d'occupation de l’état initial et de 
l’état final. En résumé, l’on peut dire que l'usage 
des statistiques quantiques (ou, ce qui revient au 
même, de la superquantification) est un élément 
aussi nécessaire à la réconciliation de la Relativité 
et des Quanta que l'introduction du spin. 


III. — Particule à spin liée. 


111.1. Transformations de Fourier cova- 
riantes [8]. — Considérons d’abord le cas, prati- 
quement si important, de la particule plongée dans 
un champ donné, invariant par translation du genre 
temps. Prenant pour repère le système propre du 
champ, il est classique de développer la solution 
générale L(x) sur le système des fonctions propres 
de l’énergie totale 


/ / D. 
W — : Æ RE : m) (38) 
2TLC LC 
sous la forme 
iprg(x)= | dés(k)p(%, k)etten, | 
è (39) 


. à 
tre (ati) 


où nous avons introduit le facteur V>r en vue de 
la suite. 


Contrairement à sa première apparence, cette 
formule est covariante relativiste. — Introduisant, 
en effet, la transformation de ‘Fourier 


(2x0 &)= ff 


= #9 


ÉEETRES 


où p désigne l'indice de spin, essentiellement non 
sommé quand il est répété, et posant 


2 


RENTE \ 
dé£e'* 1 &), | 


(40) 


= 
LC) = e(Hyn CE ki), (4) 
ainsi que 
iù=dik, idw=dir, k=h) 
x "E (4) 
(X= 199: SERA \ 
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N° 42. 
(39) devient (5) | 
i=y(e)= f]] eh et (Je) Dr. (43) 


L'intégrale de Fourier covariante (43) n’a de 
sens que si l’on définit sa norme, ce qui, en Méca-. 
nique ondulatoire, doit être fait dans l’espace. 
abstrait des fonctions de carré sommable. Postulons. 
donc, toujours sans sommation sur p, une valeur. 
finie (non invariante relativiste) pour chaque inté-. 


grale d 

Le ff ctoos= ff] vru,6e (A) 
(43) s’inverse alors suivant 

AT CE) ff] e—ihiaæ; dt (xr)du. (4) 


Ici, et dans ce qui suit, nous devrons supposer essen- 
tiellement diagonale (et hermitienne) la matrice B 
telle que 


Ÿ = d+B, T=trp 
(elle l’est en fait dans tous les cas usuels); alors, des 
valeurs finies pour les (44) entraînent une valeur 


finie pour la norme invariante relativiste 


DÉETIES 


Remarquons en passant que les seconds membres 
de (43) et (46) ont précisément la forme qui, dans 
le cas limite où le champ imposé s’évanouit, vient 
en coïncidence avec ceux des (7,) et (24), précé- 
demment données pour la particule libre. 

Maintenant se présente une difficulté heureusement 
aisée à surmonter. Physiquement, la norme n’est pas 
l'intégrale quadruple (46), mais (normalisation en 
nombre de particules) l'intégrale triple conser- 
vative, prise sur une & quelconque du genre espace 


a 
if] daidôu) = n, 
& 


formule qui, en général, n’admet pas de transformée 
de Fourier covariante. Or, il est aisé de voir qu’une 
valeur finie de la norme physique (47) entraîne 
une valeur infinie de la norme mathématique (46), 
ce qui semble renverser toute notre argumentation: 

Mais introduisons arbitrairement quatre & 
ordonnées dans le temps, &, et &, arbitrairement 
loin dans le passé, &, et €, arbitrairement loin dans 
le futur, et supposons qu'entre &, et &, d’une part, 
&; et 6, de l’autre, le champ imposé admette un 
terme additif imaginaire pur, non astreint à vérifier 
les équations différentielles du champ; par exemple, 
dans le cas d’un /{-potentiel électromagnétique, 
les A’, seront imaginaires et À, réelle, et le 4-vec- 
teur iA, n'aura pas à vérifier les équations de 


(46) 


(47) 


(5) Comparer avec Maurice Lévy [14, p. 149] qui, partant 
de (43), se ramène à une forme équivalente à (39). 


È Û pe A0 | 


— Maxwell-Lorentz. De ce fait, le 4-courant Va ne 
“sera plus conservatif entre &, et &, ni entre & 
et &,, et nous supposerons que le champ imaginaire 


est précisément tel que l'intégrale (47) tombe à 


zéro avant 6, et après &,. De la sorte, une valeur 


finie, entre &, et &;, de la norme physique (47) 


est compatible avec une valeur finie de la norme 


. mathématique (46). 


111.2. Définition covariante de l’orthogona- 


… lité. — Deux fonctions d ou € seront dites ortho- 


LE 


CRC SES A cr ae nt à ét hs Pt à + 


gonales au nouveau sens covariant si 


A da Ve do = ff} cts = 0. 


Je dis que le système des fonctions propres de 
l'énergie totale, qui est orthogonal et complet au sens 
classique, est aussi orthogonal et complet au nouveau 
sens covariant (48). 

Remarquons d’abord que, dans le 4-espace k, 
chaque fonction propre de l’énergie € (k) n’est non 
nulle que dans une bande d’énergie k,, k, + Ak, 


(A8) 


. (cas du spectre continu) ou que sur une lame infi- 


niment mince d'énergie X, (cas du spectre ponctuel), 
Soit dit en passant, ceci permet, dans le cas présent, 
d’inverser de manière covariante la formule conser- 
vative (47). Dans le 4-espace x, les transformées Ÿ (x) 
de ces € (k) ont une intensité invariante par trans- 


_ lation x,; d’une manière imagée, l’on peut dire 


qu’elles ne sont non nulles que dans des tubes indé- 


* finis. du genre temps. Chacune des précédentes 


cellules d'extension en phase, formée d’une bande 
horizontale du /-espace À et d’un tube vertical 
du 4-espace x, est occupé par n particules; mais, 
de même qu’au paragraphe II.6, il serait {out à 
fait erroné de dire qu’il y a à l'instant { n particules 
dans chaque cellule. Chaque particule franchit 


. l'instant { du temps macroscopique et c’est tout, 


sans qu’on puisse aucunement lui attribuer, dans 
cet acte, l'instant {, qui est absolument inobser- 
vable (du fait que nous parlons des valeurs de 
l'énergie). Les quatre relations d'incertitude sont 
impliquées à la fois dans le formalisme des inté- 
grales de Fourier (43), (45) et (46); la normali- 


sation physique liée à la formule (47) ne signifie 


aucunement que les particules peuvent être saisies 
sur l’hypersurface &, mais simplement qu’elles 
franchissent cette hypersurface du genre espace, 
quelle qu’elle soit. 

Ceci étant, il est clair que deux fonctions propres 
de l'énergie liées à des valeurs différentes de W 
sont orthogonales au nouveau sens (48) (du fait 
du facteur exponentiel en iAw.t). Soient donc deux € 
liés à une même valeur propre multiple W du 
spectre ponctuel. Puisque B est supposée diagonale, 
deux composantes différentes €, et €, d’un #& 
quelconque. sont orthogonales au nouveau sens 
covariant. Soient donc deux € (k) différents liés 
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Hpar 
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au même système de valeurs W et p : ils seront ortho- 

gonaux au nouveau sens suivant (49,) qui, puisque B 

est diagonale, équivaut au sens traditionnel (49,) 
GPO: EU ED: 


LT Te à 2 ge 
] (@ Co Ôt = 0; Î Gi Cp —= 0 
AT ki —E CR KQ 


(49) 


(nous avons supposé que X,; n’est pas un point 
d’accumulation du spectre ponctuel). Dans le cas 
particulier du champ de force central, les trans- 
formés Ÿ (x) des précédents € (k) sont bien connus : 
ils sont des produits d’un même facteur radial f (r) 
différentes fonctions de Laplace Y/, (6,c), 
qui forment bien un système orthogonal complet. 

Je dis maintenant que, dans le 4-espace k, les 
précédentes bandes et lames d’énergie sont symé- 
triques par rapport à l’origine des k. En effet, 
l’électron lié, par exemple, est équivalent à un positon 
lié à énergie cinétique négative, donc à une particule 
à énergie totale négative. Ceci se voit aussi en se 
reportant à la formule du spectre de l’atome hydro- 
génoide, qui comporte un radical; en considérant 
les deux signes du radical, on retombe sur ce qui 
vient d’être dit. Mais, naturellement, il serait tout 
à fait erroné d'imaginer des transitions entre niveaux 
d'énergie de signes contraires; deux niveaux symé- 
triques sont deux représentations possibles d’un 
même état physique; dans la normalisation en 
nombre de particules, c’est seulement la somme des 
nombres d’occupation correspondants qui a un 
sens. Comme pour la particule libre, l’énergie néga- 
tive est une fiction formelle suggérée par le forma- 
lisme quantique. Macroscopiquement parlant, il n’y 
a que des énergies positives. 


111.3. Passage à la superquantification. 
Représentation de Heisenberg et représenta- 
tion d'interaction. Cas d’un champ imposé 


, variable quelconque. — Tout comme ‘la for- 


mule (7) du paragraphe II, la formule (43) ci-dessus 
peut être portée directement dans le formalisme 
de Feynman [12]; on obtient ainsi la matrice de 
l’'émission-absorption dipolaire, de l'effet photo- 
électrique direct ou inverse, etc. Ici, le champ 
d'interaction induisant des transitions entre états 
libres ou liés de l’électron, entre états libres du 
photon, est le champ d’interaction entre l’électron 
et le photon libre. 

Schwinger [16], en théorie superquantifiée, a 
tiré au clair les rapports entre la représentation 
covariante de Heisenberg, et la nouvelle représen- 
tation covariante d'interaction. La représentation 
de Heisenberg est une représentation sans évolution, 
où la fonction de répartition est invariable, et où 
les équations d’onde des particules sont les équa- 
tions avec champ imposé; ceci exige un formalisme 
d’intégrales quadruples, comme on l’a vu au para- 
graphe III.2 précédent, et, par voie de conséquence, 
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à nouveau l'exclusion de toute idée d’évolution [1]. 
Au contraire, la représentation d’interaction est une 
représentation avec évolution, où le champ ne figure 
plus dans les équations d’onde des particules, mais 
induit des transitions entre états séparés de ces 
particules qu’elle met en rapport. 

Prenons un exemple. Soit un atome initialement 
et finalement libre, mais qui, entre temps, subit 
l’action d’un champ électromagnétique extérieur 
variable. Si l’on traite l’action de ce champ en repré- 
sentation d'interaction, il y a, entre l’état initial 
et l’état final, variation des poids statistiques des 
diverses valeurs des nombres d’occupation des 
états stationnaires de l’atome. Ceci est une descrip- 
tion avec évolution [1]. 

Si, au contraire, on introduit le champ extérieur 


‘variable, à côté du champ du noyau, dans l’équation 


d’onde de l’électron, ce qui (même si l’on n’introduit 
pas explicitement la fonction de répartition inva- 
riable) revient à se placer en représentation de 
Heisenberg, alors on aura affaire à une générali- 
sation simultanée des effets Stark et Zeeman (par la 
présence simultanée des deux champs, électrique 
et magnétique, et par leur variation temporelle); 
en parlant ainsi, nous avons supposé que le champ 
externe reste assez petit pour pouvoir être considéré 
comme une perturbation du champ du noyau, 
faute de quoi l’on ne pourrait plus parler de niveaux 
d'énergie, et il s'agirait d’un tout autre problème, 
dont nous dirons un mot dans un instant. Ceci est 
une description spatialisée sans évolution, impliquant 
symétriquement les quatre relations d'incertitude 
par le jeu des formules du paragraphe III.1. 
Maintenant, nous pouvons nous affranchir de 
la restriction adoptée au début du paragraphe III. 1. 
Appelons élat lié au sens large toute solution de 
équation d'onde avec champ imposé quelconque 
(arbitrairement variable). Tout le formalisme d’inté- 
grales covariantes quadruples du paragraphe III.1 
s’applique à ce cas général, y compris notamment 
la définition (48) de l’orthogonalité. Remarquons 
que l’opérateur d’ondes est self-adjoint au sens impliqué 
par ces définitions, car on a, par exemple dans le 


cas avec 4{-potentiel A}, et compte tenu du terme 
additif imaginaire que nous avons introduit au 
paragraphe III.1, 

(af o+ AE EAN) + hiÿ = 


50 
gra PAIN An — 10: (ce 


d’où, l'intégrale étant étendue à tout l'Univers, 


NI (oi 2 A7) ba ya 
= ff. . 2 ff. va 
+2 ff Ava pu; (59) 
mais (cf. $ III.1) les trois dernières intégrales sont 


simultanément nulles. Récrivant alors les (50) 
sous la forme 


|DY>= 4, AoŸ=<vD|, (2) 
il résulte de (51) que 
<y|DYy>=<vDIv), (53) 
avec par définition 
CAD = fa du. _G#) 
C, 000 


Le terme réel de masse propre k, peut être consi- 
déré ‘comme une valeur propre multiple de l’opé- 
rateur D, self-adjoint (au sens précédent), les solu- 
tions orthogonales (au nouveau sens) de l'équation 
d'onde étant les fonctions propres correspondantes. 

Ceci conduit à la description minkowskienne spatia- 
lisée d’une particule soumise à toute l'évolution 
spatio-temporelle d’un champ donné. Si l’on veut 
revenir à une description avec évolution, il suffit 
de transformer tout ou partie du champ imposé 
en terme d'interaction, et de se ramener, avec 
Schwinger [16], totalement ou partiellement, en 
représentation d’interaction. 


Manuscrit reçu le 19 mai 1954, 
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CALCUL DE LA RÉPARTITION DE L'ÉNERGIE SOLAIRE 
RÉFLÉCHIE PAR UN MIROIR PARABOLIQUE 


Par François CABANNES et ALBERT LE PHAT VINH. 


Sommaire. — L'objet de ce travail est le calcul de la répartition de l’énergie réfléchie par un miroir 
parabolique avec une source circulaire placée à l’infini ( le Soleil), dont le centre est sur l’axe du miroir; 
le diamètre angulaire est de 32° et la brillance est supposée uniforme. On détermine l’éclairement en 
chaque point d’un plan perpendiculaire à l’axe du miroir en supposant que la tache lumineuse produite 
par un point du miroir est une ellipse d’éclairement uniforme, et en intégrant les éclairements de ces 


ellipses élémentaires pour tous les points du miroir. 


Au laboratoire de Mont-Louis, l'énergie solaire 
est reçue sur un miroir concave de très grandes 


» dimensions qui est composé de nombreuses facettes 


dont la courbure et l'orientation ont été réglées 
de façon à concentrer au mieux l’énergie. La bonne 
réalisation d’un tel miroir ne peut être contrôlée 
et évaluée que par rapport au miroir parabolique 
de même distance focale et de même ouverture. 
Nous calculerons donc la répartition de l'énergie 
solaire réfléchie par un miroir parabolique au voisi- 
nage de son foyer. 


Calcul de l’éclairement en chaque point du 
plan focal. — Nous supposons que la brillance du 
Soleil est uniforme, et que son centre se trouve sur 
l’axe du miroir (ce qui est réalisé lorsque le miroir 
est en service). 

L'ouverture des miroirs utilisés est supérieure 
au double 2 f de la distance focale : dans ce cas les 
aberrations d’un objet ponctuel même voisin de 
l’axe ne sont pas calculables facilement, nous allons 
donc considérer le disque solaire dans son ensemble. 

L'énergie réfléchie par un élément de surface dS 
du miroir est 


dW = kdS cosz, 


k étant la constante solaire exprimée en watts par 
centimètre carré (compte tenu de l'absorption 
atmosphérique et du facteur de réflexion du miroir), 
et i l’angle d'incidence du rayon moyen. 

Sur le plan focal (P,) nous avons comme image 
élémentaire une ellipse (E) de centre C, dont les 
axes sont (fig. 1) : 


SE [tg (25 +16) — teg(21—16")], 


On peut supposer le centre C et le foyer F 
confondus, car pour un miroir ouvert à 2 f : 


ee Le 
FC = @igr6 -——© Z a.10o 


J=T 


(X et Y désignant les coordonnées du point S), 
Nous pouvons donc écrire 
SE ee ON Put DR EU, 
a GER 4f Fier tg 16’, 
b=(f+F)igr6. 


Fig. 1, 


Déterminons l’éclairement &; (o, Y) fourni par 
une couronne (C;) de surface AS en un point M, 
du plan focal, situé à la distance p de l’axe du 
miroir. En désignant par « l’angle solide corres- 
pondant au diamètre angulaire 32’ du Soleil, l’éclai- 
rement du point M fourni par l’élément de sur- 
face dS est (fig. 2) 
e … ÆdS cosé cosv kdS cosr 


ER ne — = ——— cos v. 
u) ie w SX2 


Deux éléments dS, et dS, symétriques par rapport 


52 
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à l’axe du miroir donnent au point M un éclaire- 
“ 6, + 6, égal au double 26, de l’éclairement du 
foyer, au tèrme en (Av) près (fi îg. 2). Pour le bord 
d’un miroir ouvert à 2f, ce terme est égal à r074X 264. 


Fig. 2. 


Si l’élément de surface dS de la couronne tourne 
autour de l’axe du miroir, l’image correspondante (E) 
tourne autour du foyer dans le plan focal et la 
surface de la couronne S,S, qui éclaire le point M 
est 
6, us Os 

27% 


DO IAIS 


) 


1 +% étant l’arc (du cercle de centre F et de 
rayon p) intérieur à l’ellipse (E) (ig. 3). 

Comme le centre de (E) est confondu avec le 
foyer, les arcs 0, et 0, sont symétriques par rapport 
à l’axe du miroir, et il est légitime d’écrire que 
l’éclairement fourni au point M par la couronne (C;) 


Fig, 3. 


est égal au produit de la surface utile de la couronne 


+ 0, 
5 = AS + 0 
par l’éclairement moyen 
kÆ cost 
& Æ 
ge rab 


d'une ellipse (E) correspondant à un élément de 
surface de la couronne dS = 1. 


# ‘CESR 
L'éclairement fourni par la couronne (Cy) est. 
donc 


&y = 


kAS cost Di+ 6 
z ab Dr 


SL dépend de la couronne (Cÿ) et 0, +0, 


dépend de cette couronne et de la distance FM = p. 
Naturellement 


Bi+ 0 920—=927x sp Z6, 


Hi+b=o sip> a. 
En intégrant &y pour toutes les couronnes du M 
miroir, nous obtiendrons l’éclairement | 


&(e) =ù Ge 


 $ 


en chaque point du plan focal (P,). Nous avons fait « 
cette intégration en partageant le miroir en 30 cou- 
ronnes. Pour chaque couronne, nous avons tracé 
AS cost 

ab 
quement 6 pour différentes valeurs de p. 

Suivant le nombre de couronnes considéré, nous 
déterminons l’éclairement & (c) pour différentes 
valeurs de l’ouverture du miroir. 


l'ellipse (E), calculé et déterminé graphi- - 


Calcul, en chaque point de l’axe du miroir, 
de l’éclairement d’une surface perpendiculaire 
à cet axe. — Les différents éléments de surface dS 
d'une même couronne donnent en un point M de 
l’axe du miroir un même éclairement, nous obtenons - 
donc sans intégration l’éclairement fourni par la 
couronne (C5) 


ÆAS cost cosv ; 4 


ave ; ; 
(0) V2 È 

= Lorxax dem PA = ? 
[(m— FR +4fFT k 

14 Aloe Fausse 4 


[(m— FR+4fYT 


m étant la distance du point M au sommet du miroir 
et l’éclairement total au point M se réduit à une 
intégrale simple 


LS ES f 
Ÿ [0] 
0 


Nous avons calculé graphiquement cette intégrale 
et nous en donnons les résultats avec ceux du calcul 
général (fig. 6). 


Y max 


CAT" 


m— 


———————— di 
[Cm — FF} + af TT 


Calcul de l’éclairement en chaque point d'un 
plan parallèle au plan focal. —— Considérons le 
plan (Py) d’ordonnées f + y. Pour une même cou- 
ronne l’image élémentaire (E') dans le plan (Py}) 


s 
É 
>." 


* 


pe. 


ke déduit de l’image (E) par une homothétie de 


centre S et de rapport fi (fig. 4). Nous 


(P,) (Py) 
Fig. 4: 


un changement d’échelle. 


F 

j 

; 

4 

1 

À pouvons donc utiliser les mêmes ellipses en faisant 
D'autre part, 


. ce qui permet de tracer le point {. 0, + 6,.est alors 
l’arc (du cercle de centre et de rayon p) inté- 
rieur à (E’). Comme les arcs 0, et 0, ne sont plus 
» égaux, nous devons supposer que l’éclairement de 
 l’ellipse (E’) est uniforme, c’est-à-dire que dans 
l'expression de l’éclairement 


TA Stoner 
EE ee cos? v, 
w S2"? 


- nous négligeons les termes en (Av). 
Cette approximation est encore légitime, car 
» pour le bord d’un miroir ouvert à 2 f ces termes sont 


. inférieurs à 2 pour 100 de l’éclairement moyen 


,  kdS cosx 
En EE < ra b. > 
L'ensemble des résultats est donné par la figure 6, 
_ qui représente en chaque point d’un plan méridien, 
 l’éclairement d’une surface plane perpendiculaire à 
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l'axe du miroir, c’est-à-dire le flux de l'énergie 


solaire réfléchie qui traverse l’unité de surface d’un 


plan parallèle au plan focal. 


Echelle. ___ 1cm 


Foyer Sens du 
qe 


axe du miroir 
parabolique 


Miroir de 6cm 
de distance focale 


ouvertures: 


La figure 7 représente pour différentes ouvertures 
du miroir la répartition de l’éclairement du plan 
focal en fonction de la distance au foyer. 
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La figure 8 représente l’énergie solaire réfléchie 
à l’intérieur d’un cercle de surface S, centré sur le 


foyer, s étant la surface de l’image de Gauss, c’est-à- 
dire 


cite de 2 bob nr 


Sr (te TON) 


Remarquons, par exemple, que si pour une ouver- * 
ture du miroir égale à o,2f l'énergie totale est « 
réfléchie à l’intérieur de l’image de Gauss, lorsque 
l’ouverture est 2 f, ce n’est plus que 60 pour 100 
de l’énergie totale qui est réfléchie sur la même 
surface. À 


Ces calculs ont été faits au Laboratoire de l'Énergie 
PNEU RS En solaire, à Mont-Louis, sous la direction de MM. Félix 
Trombe et Marc Foëx. 


1 
Dh -te-+--+--+- 


Fig. 8. : 
Manuscrit reçu le 29 mars 1954. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM. TOME 15, DÉCEMBRE 1954, PAGE 820. 


REMARQUE SUR LA 
STABILISATION DE LA RÉPONSE D'UN MICROPHONE CONDENSATEUR 


Par.C. COLIN, 


Laboratoire de Physique-Enseignement, Sorbonne. 


Sommaire, — Il est possible d'améliorer la réponse d’un microphone condensateur par contre- 
réaction si : d 


a. on adjoint une deuxième électrode; 
b. la tension recueillie est en phase avec le déplacement de la membrane. 


Nous examinerons successivement le mode de couplage utilisé et l’application qui en résulte. 


1. Notations. — a, rayon de la membrane; Point au potentiel de référence pour 


Jes potentiels variables 


o, masse spécifique superficiellé de la membrane; À 
T, tension par unité de longueur; £(r, 0), ampli- a 
tude d’un point (r, 0) de la membrane; d, distance il 
de la membrane au repos à chacune des électrodes (1) 
et (2) symétriques par rapport au plan de la mem- 
brane; 

æ(r, 0, t)= ee f), 
E,, potentiel continu de la membrane; Æ ({) poten- 
tiel de l’électrode (1); « (f) potentiel de l’électrode (2), 
dû à l'insertion d’une impédance Z, entre (2) et 
le potentiel de référence; Q,, charge de l’une 
quelconque des faces de la membrane lorsque les 
électrodes (1) et (2) sont au potentiel de réfé- 
rence; C, la capacité correspondante; Q,, Ci, Q Co Fig. 1. -- Schéma de principe. 


_-> Point au potentiel 
deréférence pour les 
potentiels continus 


-les charges et les capacités respectives des faces 
de la membrane en regard des électrodes (1) et (2); 
nous poserons 


Qi = Qo+ g1, Q:= Qo+ 9» 


et pour une couronne de surface 27r dr 


AQ@=2xrdr(5o+ ci), AQ= 277 dr(50+ 56); 


Co. 2rrdr, 


Co 27r dr 
NOTE 


Tr a? 


AC — ] NC sd 


DT 


EN NX VUS 02 
CP HNIE 
infiniment petits du premier ordre. Nous négli- 
gerons les infiniment petits du deuxième ordre. 


nous considérerons x comme des 


2. Équations de couplage. — A. ÉQUATION 
TRADUISANT L'ÉQUILIBRE MÉCANIQUE DU SYSTÈME, 
— Le déplacement Ë(r, 0,f) d’un point (r, 0) de 
la membrane tendue et maintenue fixe le long d’une 
circonférence de rayon a peut s'exprimer par une 
série double uniformément convergente. Soit 


Em09= Xi 1820 cosçn0) 


TEE EE 

+ €), (€) sin (7 0)] Yn (An r) 
avec 
Ar) 0: 


Sh(x) désigne la fonction de Bessel de première 
espèceret d'ordre n; Kmn A —Emn Son mme zéro. 
Dans le cas où l’on a uniformément 


DEEE) 0e , 
00 ris 


une simplification évidente se présente et 


AC > Tntt) To (Em). 
m>1 
A partir de cette expression, nous pouvons calculer : 


La force vive 


fi . . 
=. — 2 2 aCE ; 
d LAURE > En CE) SC Em); 


m>i 


L'énergie potentielle 
PR : V'es 
W = 2 Tr a? 2 Em Ton (Emo ); 


> m>1 


La fonction dissipation de Kelvin 


I 2 
— rs koz a? DX70 JF (Emo) 


m>i 
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k étant un coefficient spécifique pouvant être 
complexe qu'il est possible d'étudier en soumettant 
la membrane à une impulsion. 

La fonction des forces extérieures 


Dar Dino [| F(r) So(kmnr)r dr. 


mit 


Dans le cas présent nous supposons la membrane 
portée à un potentiel continu Æ;, l'électrode (1) 
à un potentiel E ({), donc au deuxième ordre près : 


F(r)= PATES +. 


Dès lors nous pouvons appliquer la méthode de 
Lagrange qui nous fournit la suite d’équations 


d'æm(t) dam) 42» _ Eo[E(t)— «(e)] 


DE ERREUR 
dé? FAGTEU ü) M Mo \ ) 2m ds PC) ( ) 
avec 
2 = w? Fi l ï S 
On = On SRatde quel que soit MI, 
CAROL T EE 
On e ‘æ? 
FE 


la quantité est homogène, au point de vue 
27T > 


e di 
équations aux dimensions, à T2. 


B. ÉQUATION TRADUISANT L'ÉQUILIBRE ÉLEC- 
TRIQUE DU SYSTÈME. — La chute de tension « se 
décompose en chutes de tension partielles corres- 
pondant aux partiels du mouvement de la mem- 
brane. Posons : 


é mo 
m>1 m>t m>1 
LE (2) 
> To (E)- 
m>A 


La membrane étant découplée, nous aurons pour 
une couronne de surface 27r dr le système d’équa- 
tions : 


AQz dglÿ AQ: 
FC RENE TE 
L 7,448 
& = Zo FPS 
Soit, en remarquant que 

AO» £s 

AC: (+ +e); 

g=ar [ CT) t)r dr 
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un système équivalent : 


Es l TS M 
De giei(e) +74 TO) 3 5, LE (D), (2) 


Ga! one SE (3) 


mA 


d l 
= 2 RO, (4) 


Fr aTque: — 2 > LE 


md 
la valeur moyenne de æx(r,t) prise à l'instant { 
sur le domaine de rayon a. 


D) n’est autre que 


(on 


3. Recherche d’une solution permanente. — 
Les équations (1), (2), (3), (4 ne sont pas indé- 
pendantes; il suffit de résoudre (1), (2), (4. 

Posons 


E(t) = 6 ejut, on, (€) —= Âm e/o?, 


a@)(#) = g*l®) ejot, a(t) = = g* ejw!, 


_ Im À = _ Ga e/ol, 


Par substitution, nous obtenons le système 
suivant : 


1 Da(Eme) S 
200 = << y (5) 


(QI 
Eo[&—jowg®], 
27 d Em Ja1( Em) J 


go [1+ 7 Cow 2 | = 


Im 


Xml or = + Jkw]—= 


g“®), œms Ym, Seront dans un rapport réel si Z, 
correspond à une capacité C'. On déduit : 


E5(6 — a*) 
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avec 


I 
LE NE 7) 
AAC ) 

m>i 


mo) = WE — &?+ Jko. 


4. Application au microphone condensateur. » 
— Nous supposons qu’une pression sonore s'exerce . 
sur la membrane au travers de l’électrode (1). 
La force s' exerçant sur une couronne de sur= 
face 27rdr devient 


EoEir dr EE 


F(r) = 


DE SE +e|+osPr dr, 


l'équation (1) devient 


d? Carat) “vdærté) 
de se dé 


+ OM Lm(Ë) 
DE) L 
p d \ Emo di ( En) 


(Lbis) 


Soit en régime permanent 


HAVE) I 1 F(&— 0* Ne 
AE En d 1 (Em) ê d | 27 d? ui D l 
à 2Co (&—x*)E0  9P* 
Er o(2 | 270 d' ps + |: 
Soit en supposant que & —— ga* (g < 0) 
_ 2 Ca Ex x) 
* ed Co Co 
= ‘ (8) 
AM Fi ( 1)(&) 
Co+ Cn Fed? SEE L 
. : 
Si | 
CHOLET SES 
Co+ Cn ex d' SANCICE PES ; | 
12 P , 
& [er 4 FE E> = E, He . 
47 dà En. 
ï E i 2 : , 
où Px;=—s (pression électrostatique). 


d 


Manuscrit reçu le 11 mars 1954. 
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me 27 £ dE SM Yi Ji(Ems) [ ne — V0)? + JA vu) | s 
+ 
a*— E àN Ja (Em) . «| 
Co nt ee 8 Dar SM A ; 
a 
Soit le résultat 
À\(S 3 
* = ( ) , ( 1 ) 
1+A(Z) 
Lord RAYLEIGH. — Theorie of Sound, chap. IX, 193 et suiv. 
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dus à une capacité interne variable, 


sous certaines influences. 


1. Distribution de vitesse des photoélectrons. 
— Dans toutes les hypothèses que nous avions for- 
mulées jusqu’à présent sur le mécanisme interne 
des photopiles, nous n'avions envisagé que le côté 
statique du problème, c’est-à-dire l’état d'équilibre. 
Nous n'avions de ce fait pas fait intervenir le temps 
mis par les photoélectrons pour traverser la couche 
d'arrêt, 

Dans lhypothèse d’Einstein, il existe un seuil 
photoélectrique correspondant à l’énergie d’extrac- 
tion de l’électron : 


hy = Rvo+ EE. 


_ Cette équation montre qu'il n’y a possibilité d’ex- 
traction d’un photoélectron que pour > »,. La 
courbe représentative ÆE —f(»), où Æ représente 
l'énergie du photoélectron, est la droite de la 
figure 1. 

Mais si nous nous plaçons à l’échelle atomique, 
ce qui revient à regarder la courbe précédente à la 
loupe, nous constatons l’absence de seuil ponctuel. 


E E 


Vo v 


Fig, r. 


Il y a raccordement, dû à la répartition statistique 
des niveaux de l’électron photoélectrique. 
* Cette répartition entraîne également l'existence 
d’une bande d'énergie autour de la courbe précé- 
dente, la loi d’Einstein étant une loi de probabilité. 
En fait, l'expérience montre que pour un faisceau 
incident de fréquence » déterminée, il existe des 
photoélectrons de vitesses très variables. 


EL 

$ 

Re. 
Mass" 


EURE 
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INTERPRÉTATION ÉLECTRONIQUE DES PHÉNOMÈNES D'INERTIE DANS LES PHOTOPILES 
(EN PARTICULIER CAPACITÉ INTERNE) 


Par G: BLET, 


Docteur ès Sciences ; 
Chef du laboratoire de photométrie du C.R.S.I.M. (Marseille), 


Sommaire. — Le calcul des durées de transit des photoélectrons à travers la couche d’arrêt permet 
d'expliquer qualitativement les phénomènes d'inertie des photopiles interprétés jusqu’à présent comme 


L'hypothèse d’une conductivité de la couche d'arrêt variable avec le nombre d’électrons libres 
donne ensuite une explication satisfaisante de la résistance interne de la photopile et de ses variations 


À partir de cet ensemble, il est possible de déterminer une constante de temps équivalente pour la 
photopile et d’en prévoir les variations en fonction des divers paramètres d’utilisation. 


La courbe de répartition de ces vitesses constitue 
la loi de distribution des vitesses scalaires. La loi 
de distribution spatiale est la loi de Lambert et la 
résultante des deux constitue la loi de distribution 
vectorielle. 

La loi scalaire que nous avons adoptée jusqu’à 


présent est une loi expérimentale conforme à la 
figure 2. 

Le maximum de cette courbe correspond à l’énergie 
la plus probable, La forme de la courbe aux points A 
et B est très mal définie. En effet, il est très difficile, 
expérimentalement, de déterminer les distributions 
d’énergie lorsque celle-ci tend vers zéro ou augmente 
indéfiniment. En particulier, on ne peut expérimen- 
talement déterminer la tangente en À ou en B, ni 
même savoir si le point B existe, c’est-à-dire si la 
courbe a un point d'arrêt ou est asymptotique 
à l’axe des V. La différence est d’ailleurs purement 
formelle et ne joue que sur le choix de l'équation 
à utiliser. 

Nous avons étudié dans un travail à part (1 les 
répercussions possibles du choix d’une formule 
convenable pour cette loi de distribution. 

Du point de vue qualitatif, seule nous intéressera 
ici l'existence de cette courbe de distribution d’énergie 
et son allure générale. 


() Voir Note C.R.S.I.M., n° 359. 
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2. Durée de transit des photoélectrons. — 
Les divers photoélectrons émis avec des vitesses 
et des directions différentes ne mettront pas tous le 
même temps pour traverser la couche d’arrêt et 
parvenir à l’électrode transparente. Le seuil photo- 
électrique des photopiles au sélénium étant situé 
dans le proche infrarouge, les énergies des photo- 
électrons expulsés par de la lumière visible sont 
déjà considérables et correspondent à des vitesses 


N 


Fig. 3. 


initiales de l’ordre de 1000 km/s environ. Mais à la 
traversée de la couche d’arrêt ces électrons vont 
perdre leur vitesse très rapidement par chocs sur 
les atomes et les ions qu’elle contient. 

Un calcul très approximatif va nous montrer que, 
même avec cette vitesse initiale considérable, le 
photoélectron met un temps appréciable, pour 
traverser une couche même très mince. 

On admet en général que le libre parcours moyen 
est indépendant de la vitesse et son ordre de gran- 
deur est 5.10 ‘em pour un semi-conducteur. 
Admettons que statistiquement l’électron perde la 
moitié de son énergie à chaque choc; sa vitesse est 
divisée par |? et le temps mis par l’électron pour 
parcourir n parcours moyens vaut : 

Dons 
1 = 


ce ul 
v K —1 


avec : 

l, libre parcours moyen -: 5.10 cm; 
D, , vitesse initiale æ 10° cm/s; 

K, facteur d’affaiblissement de vitesse. 


La durée de transit atteint o,o10 s (valeur corres- 
pondant aux constantes de temps des photopiles 
au sélénium) en donnant à n une valeur voisine de 85. 
Le chemin parcouru est donc de l’ordre de 5.10% X 85, 
soit 425 À. Cette valeur donne une limite très supé- 
rieure de l’épaisseur de la couche d'arrêt. En effet, 
les électrons n’ont pas tous des vitesses normales 
et, de plus, les libres parcours moyens successifs 
doivent s’ajouter vectoriellement et non pas sca- 
lairement. Etant donné les imprécisions sur la 
valeur de ce libre parcours moyen et sur le, facteur 
d’affaiblissement de vitesse à chaque choc, nous 
pensons que la zone de contact sélénium-métal, 
appelée couche d’arrêt, doit être contenue dans une 
épaisseur de l’ordre de 100 À à un facteur près que 
nous estimons égal à 3. 

En tout état de cause, nous admettrons donc que 
la durée de transit d’un électron à travers la couche 
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d’arrêt est inversement proportionnelle à la compo- 
sante normale de sa vitesse, sans plus. 

En procédant par ordre, à partir de la loi de 
distribution des énergies et de la loi de Lambert, 
nous pourrons calculer la loi de distribution des 
composantes normales et, à partir de celle-ci, la 
loi de distribution des durées de transit. Celle-ci 
est représentée par une courbe du type de la figure 3. 

Pour tous les calculs relatifs à ces diverses lois 
de distribution et à leurs équations, on voudra bien 
se reporter à la Note déjà citée (Note n° 359) où 
tous ces calculs sont détaillés. 


3. Établissement du courant photoélectrique. 
— À l'instant { —o, éclairons instantanément la 
cellule avec un éclairement ÆE que nous laissons 
subsister. Pendant chaque durée élémentaire dé, 
la couche photoélectrique émet un nombre de photo- 
électrons N df proportionnel à l’éclairement. 


A. Supposons qu'il n’y a pas accumulation d’élec- 
trons libres sur l’électrode collectrice, c’est-à-dire 
que les photoélectrons y arrivant trouvent un 
chemin parcouru quasi instantanément pour se 
recombiner avec les ions. positifs libérés dans Ja 
couche photoémissive (court-circuit extérieur, par 
exemple). 

Dans ce cas, à un instant { quelconque, le nombre 
total d'électrons arrivant pendant une durée dé est 
représenté par l'intégrale de la courbe des durées 
de transit prise entre zéro et £. Cette courbe a l'aspect 
de la figure 4. 


L’allure  pseudo-exponentielle de cette courbe 


Fig; 


met bien en évidence l'accroissement progressif 
du courant photoélectrique qui tend vers une limite. 

On peut définir par une courbe de ce type une 
pseudo-constante de temps comme étant le temps = 
au bout duquel la courbe n’est séparée de son 


asymptote que par la fraction : de la valeur limite. 


L'étude mathématique de ce problème indique la 
relation suivante : 
K 


| 
Vh 


rt 


7 étant cette pseudo-constante de temps; 


É 
ce 
= 

à 

f 


RICA ki 
X 


V, représentant la vitesse probable des photo- 
S électrons ; 
- À, un coefficient ne dépendant que de la forme de 
la loi de distribution des vitesses. 


B. Supposons qu'il y ait accumulation partielle 
des électrons libres sur l’électrode collectrice. L’exis- 
—._tence de charges libres crée un champ électrostatique 
à travers la couche d'arrêt. Ce champ est retardateur 
… pour les photoélectrons; son influence est double : 


…. 1° Les photoélectrons seront tous ralentis; 


DD 20 Les photoélectrons trop lents n’atteindront 
pas la couche collectrice. 


æ * 


Nr Es 


Cette seconde conséquence a déjà été étudiée à 
propos du potentiel retardateur (?). Nous avons vu 
… du point de vue statique que le résultat est une 

diminution du courant photoélectrique par rapport 

à la valeur calculée sans tenir compte de ce phé- 
|  nomèêne. 

Du point de vue dynamique, la première consé- 
quence va apporter un élément nouveau : la vitesse 
probable V, des électrons est diminuée. Si l’on 
admet, ce qui est approximativement correct, que 
la loi de distribution de vitesses n’est pas modifiée 
et que le coefficient K a conservé la même valeur, 
nous voyons que la constante de temps 7 se trouve 
augmentée en raison inverse. 

L'influence du champ électrostatique & limitée 
au premier ordre est de la forme suivante pour les 

. durées de transit 


PE, 


PRET 


lt = to + a&ti, 
d’où | 
TT + G6TR: 


Pour un même champ, la constante de temps est 
d'autant plus modifiée qu’elle est plus grande. 
En résumé, l'existence d’un potentiel retardateur 


électrique (ou la tension de sortie) et d'augmenter 
la constante de temps, c'est-à-dire la durée d’éta- 
blissement du courant. 

En fait, nous verrons qu'un autre phénomène 
beaucoup plus important, dû à la variation de 
conductivité de la couche d’arrêt, masque la variation 

de constante de temps due au potentiel retardateur. 


4. Conductivité de la couche d’arrêt. — Le 
nombre d'électrons libres contenus dans la couche 
. d'arrêt est essentiellement variable. Nous supposons 
- que dans l'obscurité, au repos, le nombre d'électrons 
. libres par unité de volume est n,. Lorsque la cellule 
est éclairée, les photoélectrons traversent la couche 
d'arrêt, mais nous supposons qu'ils n’interviennent 
+ pas dans la conductivité. Par contre, ceux qui ont 
- atteint l’électrode collectrice, retraversent cette 


(?) Voir J. Physique Rad., 1953, 14, 241. 
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a pour double résultat de diminuer le courant photo- . 
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couche d’arrêt et participent à la conduction. 
Soit n le nombre par unité de volume de ces électrons 
libres; nous représenterons la conduction de la 
couche d'arrêt par la fonction : 


= F(n+n). 


La seule hypothèse que nous faisons actuelle- 
ment sur cette fonction est qu’elle est continue et 
croissante. La résistance de la couche d’arrêt est 
donc de la forme : 

Fast K 
; "Tr Fun) 


Reprenons l'hypothèse du potentiel retardateur, 
soit : 
WE Neo), 


relation qui donne le nombre de photoélectrons 
atteignant l’électrode collectrice en fonction du 
potentiel », de cette électrode. Nous savons déjà 
que cette fonction g(v) est égale à un pour v =o 
et tend vers zéro lorsque » tend vers le potentiel 
limite 2,. 

D'autre part, la photopile étant en circuit ouvert, 
nous avons 

v= Anr, 
soit 
KA »: 


AG re et CP NOR 


(1) 

De ce groupe de relations nous allons pouvoir 
tirer quelques conclusions quant à la nature de la 
fonction F0 


19 Si nous supposons F — const., nous en dédui- 
sons que v et n sont proportionnels; comme 2 est 
limité, il s'ensuit que n l’est aussi, donc la fonc- 
tion g(v) doit tendre vers zéro lorsque l’éclairement, 
donc N, augmente indéfiniment. Donc vp,. 

Cette solution n’est donc pas impossible a priori, 
mais elle est contradictoire avec les variations de 
la résistance interne des photopiles avec l’éclai- 
rement. Nous l’éliminerons donc. 


20 Supposons ÆF variable avec n et, pour les 
valeurs très grandes de n, posons que Fest d’ordre x : 
F & n°. 

a Six < 1, nous trouvons D —> æ. Cette solution 
est à éliminer. 

DSL NT Tr AnOuS 
également à rejeter. 


trouvons v—o. Solution 


Il reste æ — 1, ce qui donne 


Ù + Vi = const. 


Cette constante n'est pas forcément égale à p,. 
En effet, la constante v, dépend de la fonction g(v), 
c’est-à-dire de la loi de distribution des vitesses, 
et de celles-ci en valeur absolue. 

Par contre, la constante v, à laquelle nous arri- 
vons dépend de la fonction de conductivité F. 


mn ARC 
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C’est évidemment la plus petite de ces deux limites 
qui sera atteinte. 

Supposons v, < 2. Dans ce cas, g() + o et n tend 
vers une limite n' telle que 


KA n° 
Mure LC + nm) 

Si, au contraire, D, < D, : go) est différent de 
zéro, n tend vers l'infini et v tend vers ».. 

Nous chercherons ultérieurement un moyen de 
différencier ces deux limites. 

De toute manière, il semble que le type de fonc- 
tion acceptable pour F soit une fonction d'ordre 1, 
la plus simple évidemment étant la loi de propor- 
tionnalité entre la conductivité et le nombre d’élec- 
trons libres. 

Ceci nous amènerait finalement à écrire pour la 
résistance interne : 

K 
ARE Ng(v)’ 
où N est proportionnel à l’éclairement. 

Ces conclusions n’impliquent pas forcément que 
les électrons libres excédentaires (n) participent 
à la conduction, maïs leur présence, en modifiant 
la charge d’espace, peut agir sur la conduction due 
à d’autres électrons libres. Nous reviendrons ulté- 
rieurement sur ce point. 


5. Influence combinée de la durée de transit 
des photoélectrons et de la variation de conduc- 
tivité de la couche d'arrêt. — Nous avons admis 
pour la courbe d’arrivée des photoélectrons une 
courbe de forme exponentielle à constante de 
temps définie. Prenons cette constante de temps 
pour unité et écrivons : 


n = N(i—et), où N—KE, 
le nombre de photoélectrons libérés étant pro- 
portionnel à l’éclairement. 

D'autre part, nous avons admis pour la résistance 
interne une expression de la forme : 


A 
ji + 


A — 


en négligeant pour l'instant l'influence du potentiel 
retardateur. 
En régime transitoire, à chaque instant, la tension 
aux bornes sera | 
DNA RAI TESTÉES (2) 
no + N(1—-eTt) ; 


[en tenant compte du potentiel retardateur, il fau- 
drait remplacer N par N qg(v)]. 

Cette expression (2) exprime la croissance de la 
tension aux bornes lorsqu’à l'instant { —o est 


appliqué instantanément un éclairement E. L’expres-. 
sion peut aussi se mettre sous la forme : 


FC) 


où / (#) est la fonction caractéristique d’arrivée des L; 
photoélectrons (1 — e-!, par exemple) et . le rapport f 
du nombre d'électrons libres à l'obscurité au nombre # 
d'électrons libres excédentaires par unité de temps, « 


dus à l’éclairement Æ: 
En régime permanent, f({) tend vers l’unité et 


RUES. 3 : N 
Donc pour les éclairements très faibles : v w D, = 4 
‘ 0. 


est proportionnel à l’éclairement. 
Au contraire, pour les éclairements très élevés, 


4 


Fis,n5: 


le pôtentiel croît lentement vers la valeur asympto- 
tique D. 

Pour l'instant, nous nous intéressons particuliè- 
rement au régime transitoire et il est aisé de voir 
que les diverses courbes de croissance de » sont 
différentes et d’autant plus rapides que l’éclairement 
est plus élevé. | 

Nous avons tracé sur le graphique de la figure 5 
les courbes de croissance de la tension pour 


UN) 


N 


—=10—1<—0,1—0,01, 


en ramenant chaque fois la valeur asymptotique à 
la même ordonnée pour comparer plus aisément 
les formes de ces courbes. 

Lorsque l’éclairement tend vers zéro, la forme de la 
courbe tend vers celle de f (#). 

Sur cette figure, pour simplifier les calculs, nous 
avons pris 

J(t) =1— €", 


ce qui donne par définition la valeur 1 pour la 


constante de temps aux éclairements très faibles. 
En prenant la définition classique de la pseudo- 
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constante de temps, nous trouvons pour ces diverses | 
- courbes les constantes de temps suivantes : 


= CO (ce) 10 I 0,1 O,01 


0,94 0,618 0,017 


Cette variation considérable visible sur le gra- 


ES phique est bien mise en évidence expérimentalement. 
\ Les photographies ci-jointes le prouvent qualita- 
. tivement (nous poursuivons la mise au point d’un 


montage devant nous permettre des vérifications 


quantitatives). 


D'autre part, on peut se rendre compte de la 
variation du niveau de la tension de sortie en fonc- 
tion de l’éclairement. Cette variation telle qu’elle 
est déduite de l’équation établie ci-dessus peut se 
traduire par une expression de la forme : 


forme déjà rencontrée et ne représentant correc- 
tement les variations de potentiel que pour les 
éclairements faibles. En effet, nous avons négligé 
l'influence du potentiel retardateur qui devient 
plus importante aux éclairements élevés. 

Dans ce cas la formule initiale devient : 


ft) 


+ f(t) 


DV = Vo 
70 


N g(v) 
d’où l’on tire une expression de la forme 


E = A LE PLASE EPS 0 

(vo—v)g(v) 

Cette expression ressemble beaucoup à celle que 
nous avions admise jusqu’à présent en tenant compte 
seulement du potentiel retardateur : 


Ù 
FE = A — . 
g(v) 
r. ” ch I ô 
Le terme correctif D VOISIN de 7 aux faibles 
AN 0 


éclairements, ne modifie donc pas cette région de 
la courbe. Pour les éclairements élevés, par contre, 
ce terme augmente la valeur de l’éclairement néces- 
saire pour avoir une tension donnée. 

Dans le domaine des éclairements réalisables et 
mesurables au laboratoire, les deux courbes sont 
très voisines l’une de l’autre : en effet, la cons- 
tante À ne peut se déterminer avec précision que 
pour les éclairements très faibles, et, dans ce cas, 
la tension de sortie devient elle-même trop faible 
pour être mesurée avec précision. 

La formule (2) semble donc représenter correc- 
tement à la fois le phénomène transitoire et le phé- 
nomène permanent. 
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Voyons toutefois quelle est l'influence du fac- 
teur g (v) sur le transitoire. 

Le début de la courbe est très peu modifié puisque v 
est très petit, donc g(v) très près de l'unité. Par 
contre, pour l'extrémité de la courbe, g(v) diminue 
et tend vers une valeur g(v;) plus petite que l'unité, 


Fig. 6. 


donc la limite atteinte v, est plus petite que celle 
qui serait atteinte sans introduire g (v). 
Le calcul de la nouvelle pseudo-constante de 
temps donne les résultats suivants : 
uote Le) 10 I 


0,1 0,01 


0,80 0,445 0,179 0,054 
d’où il résulte que l'effet du potentiel retardateur 
se traduit par une diminution de la constante de 
temps pour les faibles éclairements et une augmen- 
tation d’autant plus marquée que l’éclairement 
est plus intense, pour les forts éclairements. Le 
phénomène prépondérant reste toutefois la dimi- 
nution de rt avec l’éclairement. 
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6. Influence de la longueur d'onde. — Nous 
avons vu précédemment que deux possibilités 
s’offraient à nous pour l'interprétation du potentiel 
limite. L'une des valeurs était conditionnée par la 
loi de distribution de vitesses des photoélectrons, 
l’autre par la loi de conductivité de la couche 
d'arrêt. La première de ces deux valeurs peut 
dépendre de À alors que la seconde n’en dépend 
certainement pas. 

Or, expérimentalement, il ne nous a pas été pos- 
sible de mettre en évidence, aux erreurs d’expé- 
rience près, une influence quelconque de la longueur 
d'onde de la lumière excitatrice sur le potentiel 
limite. Il semble donc que ce soit la fonction de 
conductivité de la couche d’arrêt qui régisse « ce qui 
se passe » dans la photopile. 

Une autre confirmation nous est apportée par le 
fait que la longueur d’onde est également sans 
influence sur la constante de temps à condition de 
choisir pour chaque À un éclairement produisant 
le même courant de court-circuit. Les phénomènes 
sont donc conditionnés uniquement par le nombre 
d'électrons libérés, indépendamment de leur énergie 
initiale. 


Conclusions. — En tenant compte d’une part 
de la loi de distribution de vitesse des photoélectrons 
d'autre part de la variation de conductivité de la 
couche d’arrêt en fonction du nombre d'électrons 
libres la parcourant, enfin, de l’effet du potentiel 
retardateur, il semble possible de représenter l’en- 


CR RS RE SRE mine 1 0, 7 
(à ‘ ÿ À Le We Cv" Cru or 
FAUrE 


N°12. 
semble des phénomènes transitoires et permanents | 
dans une photopile par une relation de la forme : 

ONE 
+ f(E) 


Ùv = Vo A 

E g(v) 

où : 

À est une constante; 

E est l’éclairement; 

g @), la fonction du potentiel retardateur du type 
vd \ 

ErT 

f(, la fonction de croissance du courant d’élec- 


trons libres, peu différente en pratique 
L 
detre peus 
v le potentiel limite caractéristique du type 


de couche d’arrêt utilisé. 


I (n =, en pratique); 


Il semble confirmé que l'énergie initiale des 
photoélectrons n'intervient pas dans cette relation 
qui reste valable en lumière non monochromatique, 
à condition toutefois que la composition de la 
lumière ne se modifie pas lorsque son intensité varie. 

L'expérience confirme ces résultats qualitative- 
ment. Pour une vérification quantitative, nous 
attendons la mise au point d’un montage délicat 
qui doit nous permettre un relevé très précis des 
transitoires et des mesures assez précises des pseudo- 
constantes de temps. 


Manuscrit reçu le 8 mai 1954. 
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LETTRES A LA RÉDACTION 


ÉTATS DE SURFACE ET CONDENSATION 
DE JETS MOLÉCULAIRES 


Par P. GarRiN et P. PRUGNE. 


Récemment, Devienne [1], [2] signalait l’importance 
de la nature du support récepteur sur la condensation 
des jets moléculaires à la température normale. Il 
semble indiquer que l’état de surface du support 
a moins d'influence que sa nature [1]. 

Notre expérience [3], [4], [5] nous incite à penser 
que l’état de surface (pureté, dégazage, etc.) est 
d’une importance majeure et que, s’il n’est pas 
rigoureusement défini, il masque complètement l'effet 


_ dû à la nature de la surface. 


Nous avons étudié la condensation du zinc et du 
cadmium sur aluminium nu, sur aluminium recouvert 
de sa couche naturelle d’alumine, sur calcium, cuivre, 
argent, etc. 


Méthode et dispositif expérimentaux utilisés. 
— Les feuilles a, b, c (fig.) sont en aluminium recou- 
vert de son alumine, nettoyées classiquement par 
solvants (alcool pur) et pesées. Un écran d permet 
de masquer aux rayons moléculaires les feuilles a 
ou c. Le creuset évaporateur 1, l’écran d masquant c, 
permet d’obtenir sur a et b le support sur lequel on 
veut étudier la condensation de l’élément considéré. 
Les dépôts sur a et b sont identiques par raison de 
symétrie. L’écran vient alors masquer a. A la fin 
de l'expérience, la pesée de a donne la masse super- 
ficielle de ce premier dépôt. 

Le creuset 2 permet l'obtention des jets molé- 
culaires étudiés. 

Des conditions géométriques [6], nous déduisons 
les masses totales M, et M. de l’élément qui viennent 
frapper les surfaces b et c. Par pesées, nous détermi- 


_ nons les masses effectivement condensées my, m,… 


Ces mesures donnent les coefficients de conden- 
sation : 


m 2V2 12,2 z 
K = de l'élément considéré sur support pré- 
My 
évaporé sous vide; 
ST de l'élément considéré sur Al recouvert 
[a 


de sa couche d’alumine. 


Nous trouvons : 


Zn sur Al pré-évaporé K — o sur Al recouvert d’alu- 
mine k = 1; 

Zn sur Ca pré-évaporé K = 1,6 sur Al recouvert 
d’alumine k = 0,6; 

Zn sur Cu pré-évaporé K + 1,5 sur Al recouvert 
d’alumine k = 0,5; 
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Zn sur Ag pré-évaporé K = 0,9 sur Al recouvert 
d’alumine k = 0,75; 

Cd sur Al pré-évaporé K = 1 sur Al recouvert d’alu- 
mine X = o. 


Conclusions. — Les résultats trouvés n’indi- 
quent que les valeurs relatives des coefficients de 
condensation. Néanmoins, on peut déjà en tirer les 
conclusions suivantes : 


19 Les facteurs de condensation pour le zinc et 
le cadmium sont très caractéristiques de l’importance 
de l’état de surface. La valeur K = : pour le cad- 
mium sur aluminium était en particulier passée 


Différents Supports 
de condensation 


Ecran évitant 
Ja vue directe 
entre les deux 
Sources 


Orifice 
de pompage 


, Creuset 
émetteur 2 


l'écran 


Coupe ab 


inaperçue [2], de même que l’espèce d’effet de répul- 
sion de l’aluminium pur sur le zinc. 

Il est admis de même que KX — o pour le cadmium 
se condensant sur verre [2] à la température normale. 
Or Cd se condense avec K = r sur du verre qui a été 
nettoyé par effluvage électrique dans l’air. 

20 Les valeurs de X sur Cu et Ca sont paradoxa- 
lement supérieures à 1. Ceci peut s’expliquer par le 
fait que les molécules émises en dehors de l’angle 
solide déterminé par le support et le creuset, après 
plusieurs réflexions sur les parois de la cloche, sont 
venues se condenser sur ces supports, de Cu et Ca. 
Ainsi, certains métaux pré-évaporés sous vide jouent 
un rôle de piège sur les faisceaux moléculaires par 
rapport à d’autres métaux. Ceci pourrait, semble-t-il, 
être utilisé pour la préparation de corps rares. 

30 Enfin, tous les dépôts sur supports pré-évaporés 
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sous vide ont un aspect bien plus uniforme et homo- 
gène que ceux obtenus jusqu’à présent sur Al, Cu 
ou verre. 


Ces travaux ont été exécutés sous la direction de 
M. Berthelot, chef du Service de Physique nucléaire. 


Manuscrit reçu le 4 juillet 1954. 
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BANDES NOUVELLES DANS LE SPECTRE D'ÉMISSION 
DE L'ACÉTALDÉHYDE 


Par Mme RÉGNIER, 
Chimie théorique, Faculté des Sciences, Nancy. 


Nous avons obtenu un spectre très complexe par 
décharge haute fréquence sans électrodes dans un 
courant gazeux d’acétaldéhyde. Les conditions d’expé- 
riences et les premiers résultats ont fait l’objet de 
deux publications [1]. Ce spectre contient surtout 
les produits de décomposition de la molécule d’acé- 
taldéhyde, mais aussi un certain nombre de bandes 
qui semblent n’avoir jamais été signalées. Ces bandes 
sont dégradées vers le rouge et s’étendent entre 3 7950 
et 5 ooo À. Les positions approximatives des prin- 
cipales têtes de bandes sont les suivantes : 4 7928, 
4715, 4 699, 4 300 (?), 4 144; 3 789, 3 782,3 767 
CLA9:708 27 

L’extrême complexité du spectre n’autorisant pas 
des recherches plus précises, nous avons abandonné la 
décharge haute fréquence trop violente et donnant 
donc trop de produits de dissociation diatomiques, 
et nous avons employé une décharge dans une cathode 
creuse. Les électrodes sont inspirées d’un montage 
de Schüler [2]. Par ailleurs, les conditions d’expé- 
riences sont sensiblement les mêmes que dans le cas 
de la décharge haute fréquence, à l’exception de la 
pression, plus élevée que précédemment. 

Des diverses observations ainsi faites, il ressort que : 

19 Les bandes inconnues, assez intenses dans la 
décharge haute fréquence, n'apparaissent absolument 
pas dans la cathode creuse; 

2° Le spectre de CO, très développé dans la décharge 
haute fréquence est peu intense dans la cathode 
creuse : seules apparaissent les bandes d’Angstrôm. 
Les bandes de triplets ont disparu totalement; 

3° Le spectre du carbone (bandes de Swan), 
‘absent dans la décharge sans électrodes, est intense 
dans le second cas; 

4° Les bandes de CH sont très fortes dans la cathode 
creuse, alors qu’elles n’apparaissent — très faible- 
ment — dans la décharge haute fréquence que lorsque 
la pression augmente un peu. 


Ainsi les bandes inconnues sont apparues en même 
temps que les bandes de CO et en l'absence des 
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bandes de C, et de CH. On est donc en présence de 
trois hypothèses possibles sur l’origine de ces bandes : 

— Elles peuvent être un nouveau système de CO 
non encore signalé; 

— Elles peuvent être dues à la molécule d’acé- 
taldéhyde elle-même (cette hypothèse est sans doute 
la moins probable, bien que récemment un spectre 
d'émission de l’acétaldéhyde lui-même ait été signalé 
dans la région 3 500-4 200 À [3]; 

— Elles peuvent avoir pour émetteur un grou- 
pement plus ou moins complexe formé à partir de 
carbone, hydrogène et oxygène. 


C’est ce problème que nous nous efforçons de 
résoudre en étudiant maintenant la structure vibra- 
tionnelle et surtout la structure fine des bandes inté- 
ressantes. 
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VARIATIONS DE LA VITESSE DES ULTRASONS 
DANS LE PROPANE EN FONCTION DE LA TEMPÉRATURE 
ET DE LA PRESSION 


Par A. LACAM, 
Laboratoire des Hautes Pressions, Bellevue, 


Utilisant le même appareillage qui nous a permis 
de faire des mesures dans le méthane sous pression [1] 
jusqu’à 200° C, nous avons mesuré la vitesse de propa- 
gation des ultrasons dans le propane pour les tempé- 
ratures de 25, 50, 96,8, 125, 150, 175, 200 et 22500 
et des pressions allant jusqu’à environ 1100 atm. 
Le propane que nous avons utilisé nous a été fourni 
par la Phillips Petroleum (U.S. A.) qui le garantit 
spectroscopiquement pur dans l’infrarouge. Les tempé- 
ratures expérimentales étaient maintenues constantes 
au 1/20 de degré. La température critique du gaz 
étant 960,8, nous avons pu opérer dans la phase 
liquide (isothermes 25 et 500 C), région dans laquelle 
le fluide a un comportement normal : la vitesse 
décroît avec la température. La région des pressions 
peu élevées, en particulier celle du point critique a été 
explorée et représentée sur la figure pour donner une 
vue d’ensemble des phénomènes. Cette région a été 
étudiée d’une façon plus systématique par J. Noury [2]. 
On note un bon accord entre les deux séries de résul- 
tats. Dans la phase gazeuse nous observons que la 
vitesse croît avec la température pour les faibles 
pressions; quand la pression augmente il se produit 
une inversion de la position relative des isothermes. 
Ainsi pour des pressions supérieures à 120 atm, la 
vitesse décroît avec la température. Aux pressions 
élevées l’ensemble des isothermes compris entre 25 
et 2250 C présente l’allure d’un faisceau de paraboles 
sensiblement parallèle. Les courbes isobares V = /(£2) 
nous ont permis de déduire l’isotherme 75°C, qui 
est représentée en traits interrompus sur la figure. 
Nous n’avons décelé aucune dispersion dans le domaine 
de fréquences étudiées, qui s’étend entre 3 et 4 MHz. 


id 


 Rowlinson [3] qui n’ont pas trouvé de dispersion 
aux faibles pressions pour des fréquences comprises 


V m/s 


1500 | 
1400 


LETTRES A LA RÉDACTION 


Ce résultat est à rapprocher à celui de Lambert et 


LT OR 


831 
entre 0,25 et 4 MHz. La précision sur les mesures 


de la vitesse est de l’ordre de 1 pour 100 (sauf dans 
la région critique et aux pressions inférieures à 80 atm). 


Manuscrit reçu le 15 octobre 1954. 
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VITESSE DES ULTRASONS DE 585 kHz 
DANS LE PROPANE AU VOISINAGE DE L’ÉTAT CRITIQUE 


Par J. Noury, 
Laboratoire des Hautes Pressions, Bellevue. 


Poursuivant mes -études concernant la détermi- 
nation de la vitesse des ultrasons dans les fluides 
autour de la région critique par la méthode de la 


600 700 800 300 1000 


Patm. 


diffraction de la lumière, j'ai effectué des mesures 
dans C3Hy | 

L'installation utilisée précédemment pour l’étude 
des régions critiques de CO,, N,0 et CH, [1], [2] a 
été considérablement modifiée. La bombe expéri- 
mentale a été placée cette fois dans un double bain 
de polyéthylène-glycol de dimensions notables for- 
tement brassé. Ceci a permis d’étendre le domaine 
des températures explorées. 

La vitesse des ultrasons de fréquence 585 kHz a 
été mesurée le long des isothermes 75, 90, 96,80 C 
(isotherme critique); 100 et 1250C maintenues 
constantes au 1/100€ de degré pour un domaine de 
pression allant de 10 à 150 atm. Le propane provenait 
de la Phillips Petroleum qui le garantit spectroscopi- 
quement pur dans l’infrarouge. 


832 


La précision des mesures est 1,5 pour 100 dans la 
région critique. Les résultats obtenus pour les pres- 
sions plus élevées sont en bon accord avec ceux de 
A. Lacam [3]. 

Comme dans les études précédentes, on constate 
que lorsque la pression croît la vitesse décroft nota- 
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blement, passe par un minimum au point critique : 
puis croît à nouveau. La vitesse minimum énre- | 
gistrée sur l’isotherme critique (96,89 C) est 152,5 m/s . 
pour 42,6 atm; c’est la valeur la plus proche au “point 
critique qui ait pu être mesurée. 

La stabilité du propane soumis aux ondes ultra- 
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sonores a été telle que des spectres de difiraction 
relativement nets ont pu être observés dans la région 
critique et un grand nombre de points ont été relevés 
dans cette région à des intervalles très rapprochés. 

En conséquence l’isotherme critique présente un 
caractère de continuité bien plus accusé que dans les 
études précédentes, 

Sur l’isotherme 75° C on constate un accroissement 


Le Gérant : MAURICE BLONDIN. 
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très brusque de la vitesse autour de la pression 30 atm, 
ce qui correspond au passage de la phase vapeur à 
la phase liquide. 
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